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ABSTRACT 


'Results  are  preseoted  from  a  towing-tank  experiment  of  propeller-hull 
interaction  conducted  in  order  to  provide  detailed  documentation  of  the 
complete  flow  field  appropriate  both  for  explicating  the  flow  physics  and 
validating  computational  methods.  Mean-velocity  and  pressure  field  meas¬ 
urements  were  made  for  the  with-  and  without-propeller  conditions  for  the 
Series  60  Cg  =  .6  hull  form  at  numerous  stations  both  upstream  and  down¬ 
stream  of  the  propeller  and  in  the  near  wake  region.  Surface-pressure 
distributions  and  wave  profiles  were  measured  for  both  conditions.  Resis¬ 
tance  and  self-propulsion  tests  were  also  conducted.  The  experimental 
equipment  and  procedures  are  described,  and  the  results  are  discussed  to 
point  out  the  essential  differences  between  the  flows  with  and  without 
propeller.  The  results  are  analyzed  to  assess  the  nature  of  the  interac¬ 
tion  between  the  propeller  and  the  hull  boundary  layer  and  wake.  To  this 
end,  use  is  made  of  a  propeller-performance  program  with  both  nominal  and 
effective  Inflows.  It  is  shown  the  most  features  of  the  interaction  can 
be  explained  as  a  direct  consequence  of  the  propeller  loading  resulting 
from  its  operation  with  a  three-dimensional  nonuniform  inflow. 
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I.  INTRODDCTION 


Detailed  experimental  information  concerning  the  interaction  between  a 
propeller  and  the  flow  over  a  ship  hull  is  very  limited,  especially  for 
practical  hull  forms.  The  most  extensive  data  available  are  those  con¬ 
cerned  with  simple  propeller-shaft  configurations.  These  include  steady 
(circumferentially-averaged)  mean-velocity  profiles  in  the  immediate  vicin¬ 
ity  of  the  propeller,  and  in  the  near  and  intermediate  wake  regions,  and 
some  limited  unsteady  (phase-averaged)  mean-velocity  profiles,  and  steady 
and  unsteady  turbulence  profiles  (see  Stern  et  al.  (1988b)  for  references). 
The  situation  for  axlsymmetric  bodies  is  similar,  although  the  data  are 
considerably  more  limited.  Hucho  (1968)  and  Huang  et  al.  (1976,  1980) 
report  wind-tunnel  measurements,  and  Nagamatsu  et  al.  (1978)  and  Toda  et 
al.  (1982)  describe  towing-tank  measurements,  of  steady  mean-velocity  pro¬ 
files  upstream  of  the  propeller.  All  these  cases,  except  Huang  et  al. 
(1976,  1980),  precluded  realistic  wake  measurements  because  the  propeller 
was  driven  from  downstream.  Schetz  and  associates  (1975,  1981,  1983)  and 
Neu  et  al.  (1988)  report  wind-tunnel  measurements  of  steady  mean-velocity 
and  turbulence  profiles  in  the  near  and  intermediate  wake,  including,  in 
some  cases,  the  effects  of  appendages  and  angle  of  attack.  Rood  and 
Anthony  (1988)  investigate  propeller-appendage  interaction  through  detailed 
steady  (i.e.  ,  averaged  over  a  sufficient  number  of  propeller  revolutions) 
and  unsteady  mean-velocity  measurements  in  a  plane  just  upstream  of  the 
propeller.  Experiments  have  also  been  performed  for  a  vertical  flat  plate 
with  a  propeller  operating  near  the  trailing  edge  driven  from  downstream  in 
which  steady  mean-velocity  profiles  were  measured  upstream  of  the  propeller 
both  in  wind  tunnels  (Hucho,  1968)  and  in  towing  tanks  (Toda,  1984; 
Nagamatsu,  1985).  The  latter  reference  includes  results  for  T  sections. 
Similar  experiments  have  also  been  performed  for  the  Wigley  parabolic  hull 
at  low  Froude  number  by  Toda  et  al.  (1984)  and  Sato  et  al.  (1986). 

Although  numerous  experiments  have  been  performed  in  towing  tanks  fcr 
various  practical  hull  forms,  none  is  sufficiently  detailed  to  document  the 
entire  flow  field.  All  are  focused  on  some  particular  aspects  of  the  gen¬ 
eral  problem  of  propeller-hull  interaction.  The  most  notable  recent  exper¬ 
iments  have  been  those  performed  in  Japan  in  support  of  the  development  of 
ultra-large  energy-efficient  merchant  ships,  i.e.,  low-speed,  large-block- 
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coefficient  tankers  and  ore  carriers,  and  high-speed  container  ships  (see, 
e.g.,  Tasaki  et  al.,  1987).  Typical  studies,  citing  example  references, 
are  as  follows.  Steady  mean-velocity  profiles  have  been  measured  at  one  or 
two  stations  upstream  of  the  propeller  for  the  with-  and  without-propeller 
conditions  to  study  scale  effects  (Dyne,  1974;  Kux  and  Laudan,  1984), 
unstable  phenomena  encountered  during  self-propusion  tests  (Taniguchi  and 
Watanabe,  1969),  and  the  characteristics  of  effective  wakes  (Laudan, 
1931;  Kasahara,  1985,  1986).  Similar  measurements  have  been  made  in  the 
near  and  intermediate  wake  to  study  recovery  of  the  rotational  energy 
(Ishida,  1986)  and  propeller-rudder  interaction  (Baba  and  Ikeda,  1980). 
Also,  steady  surface-pressure  distributions  have  been  measured  to  study 
thrust  deduction  (Cox  and  Hansen,  1977;  Fuji!  and  Fukuda,  1984).  Lastly, 
only  two  studies  report  limited  steady  turbulence  measurements  (Chen,  1964; 
Kux  and  Laudan,  1984). 

It  is  apparent  from  the  foregoing  that  experimental  studies  of  propel¬ 
ler-hull  interaction  are  required  to  provide  detailed  documentation  of  the 
complete  flow  field  appropriate  both  for  explicating  the  flow  physics  and 
for  validating  computational  methods.  The  present  cooperative  study  be¬ 
tween  The  University  of  Iowa  and  Osaka  University  was  undertaken  for  this 
purpose.  In  particular,  steady  mean-velocity  and  pressure  field  measure¬ 
ments  were  performed  for  the  with-  and  without-propeller  conditions  for  the 
Series  60  Cg  =  .6  hull  form  at  numerous  stations  both  upstream  of  the  pro¬ 
peller  and  in  the  near  wake  region.  Surface-pressure  distributions  and 
wave  profiles  were  measured  for  both  conditions.  Also,  resistance  and 
self-propulsion  tests  were  conducted.  The  experiments  were  performed  in 
the  Osaka  University  towing  tank  at  low  Froude  number,  Fr  =  .16,  to  mini¬ 
mize  free-surface  effects.  The  test  conditions  and  results  are  documented 
in  sufficient  detail  to  be  useful  as  a  test  case  for  validating  computa¬ 
tional  methods. 

Two  4m  long  models  were  constructed  for  the  experiments:  a  wooden 
model  used  for  the  mean-velocity  and  pressure  field  measurements  (figure 
la);  and  a  fiber-reinforced  plexiglass  model  with  pressure  taps  used  for 
the  surface-pressure  measurements  (figure  lb).  A  145.64mm  diameter,  5 
bladed  propeller  with  MAU  sections  was  used.  The  with-propeller  measure¬ 
ments  were  performed  for  the  model  self-propulsion  condition.  The  princi- 
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pal  dimensions  of  the  models  and  propeller  are  given  in  Table  1  and  the 
locations  for  the  detailed  measurements  are  summarized  in  Table  2. 


In  the  presentation  of  the  results  and  the  discussions  to  follow,  a 
Cartesian  coordinate  system  is  adopted  in  which  the  x-,  y-,  and  z-axes  are 
in  the  direction  of  the  uniform  flow,  port  side  of  the  hull,  and  increasing 
depth,  respectively.  The  origin  is  at  the  intersection  of  the  waterplane 
and  the  forward  perpendicular  of  the  hull.  The  mean-velocity  components  in 
the  directions  of  the  coordinate  axes  are  denoted  by  (u,v,w)  and  the 
carriage  velocity  by  U.  Unless  otherwise  indicated,  all  variables  are 
nondimensionalized  using  the  model  length  between  perpendiculars  L,  car¬ 
riage  velocity  U,  and  fluid  density  p. 

The  Series  60  Cg  =  .6  hull  form  was  selected  for  the  experiments  as  a 
representative  fine  hull  form  and  to  complement  the  many  previous  and  on¬ 
going  studies  with  this  geometry.  The  Series  60  geometry  was  conceived  to 
provide  systematic  Information  on  the  design  of  lines  for  single-screw 
merchant  ships  ca.  1950.  The  parent  form,  Cg  =  .6,  was  designed  based  on 
considerations  of  then  successful  ship  designs.  A  full  account  of  the 
original  methodical  series  is  provided  by  Todd  (1963).  The  many  experi¬ 
ments  with  the  Series  60  hull  form  since  then  are  far  too  numerous  to  re¬ 
view;  however,  of  particular  interest  is  the  fact  that  it  is  one  of  the 
four  hull  forms  selected  for  the  Cooperative  Experimental  Program  (CEP)  of 
the  Resistance  and  Flow  Committee  of  the  International  Towing  Tank  Confer¬ 
ence  (ITTC,  1987).  The  experiments  under  this  program  are,  however,  re¬ 
stricted  to  hulls  without  propellers. 

An  outline  of  this  report  is  as  follows.  The  experimental  equipment 
and  procedures  are  described  in  Chapter  II.  In  Chapter  III,  the  results 
are  presented  and  discussed  to  point  out  the  essential  features  of  the  flow 
for  both  the  with-  and  without-propeller  conditions.  Next,  in  Chapter  IV, 
the  results  are  analyzed  to  assess  the  nature  of  the  interaction  between 
the  propeller  and  the  hull  boundary  layer  and  wake.  Lastly,  in  Chapter  V, 
some  concluding  remarks  are  made. 
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II.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


The  experiments  were  performed  In  the  Osaka  University,  Department  of 
Naval  Architecture,  towing  tank  (figure  2).  The  towing  tank  is  100m  long, 
7.8m  wide,  a^'J  4.35m  deep.  The  towing  carriage  is  driven  by  four  15kw  DC 
motors  whi  ..i  are  controlled  by  a  static  Leonard  system  with  analog-digital 
hybrid  control.  The  carriage  can  be  driven  at  speeds  up  to  3.5m/s  with  an 
accuracy  of  ±  Imm/s. 

A.  Ship  Models  and  Propeller 

The  lines  of  the  Series  60  Cg  =  .6  ship  models  used  in  the  experiments 
are  shown  in  figure  3.  These  conform  to  the  standard  offsets;  however,  a 
small  modification  was  made  to  the  stern  geometry  (.955  x  ^  .995)  in 
order  to  attach  the  propeller.  The  broken  lines  in  figure  3  show  the  ori¬ 
ginal  bare  hull  and  the  solid  lines  the  modified  one.  The  details  of  the 
stern  arrangement,  including  the  modification  to  attach  the  propeller  and 
the  location  of  the  propeller  plane,  were  based  on  the  original  methodical 
series  (Todd,  1963).  The  principal  dimensions  of  the  models  are  given  in 
Table  1  and  the  offsets  are  provided  in  Table  3.  Two  4m  long  models  were 
constructed  for  the  experiments:  a  wooden  model  was  used  for  the  mean- 
velocity  and  pressure  field  measurements  (figure  la);  and  a  fiber-rein- 
f orced-plexlglass  model  with  pressure  taps  was  used  for  the  surface-pres¬ 
sure  measurements  (figure  lb).  In  order  to  induce  turbulent  flow,  a  row  of 
trapezoidal  studs  with  1.5mm  height,  1.5mm  length,  2mm  front  width  and  1mm 
back  width,  were  fitted  at  lOmm  spacing  on  both  models  at  x  =  .05. 

The  lines  of  the  propeller  used  in  the  experiments  are  shown  in  figure 
4a.  The  principal  dimensions  are  given  in  Table  1  and  the  offsets  are  pro¬ 
vided  in  Table  4.  This  is  a  conventional  stock  propeller  designed  based 
on  the  MAU  methodical  series  (Tsuchida  et  al.  ,  1958)  with  145.64mn  dia¬ 
meter,  constant  pitch,  zero  skew,  6  degree  rake,  5  blades,  and  MAU  n  =  25 
sections.  Although  there  are  some  minor  differences,  here  again,  the  di¬ 
mensions  were  selected  based  on  the  original  methodical  series  (Todd, 
1963) . 
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B.  Instrumentation 


Single-  and  duplex-balance  rod-deflection  type  dynamoneters  were  used 
for  the  resistance  and  propeller  open-water  tests,  respectively.  These 
same  dynamometers  were  also  used  for  the  self-propulsion  tests. 

The  surface-pressure  distribution  was  measured  with  359  Inm-diameter 
pressure  taps  located  in  the  stern  region  (x  _>.  where  the  influence  of 
the  propeller  was  expected  to  be  large  (figure  5a).  The  holes  were  spaced 
sufficiently  close  in  both  the  girthwise  (for  9  stations)  and  axial  (for  3 
waterlines)  directions  to  allow  for  evaluation  of  pressure  gradients.  In 
order  to  ascertain  the  effects  of  hole  interference,  the  starboard  spacing 
was  made  somewhat  coarser  than  the  port.  Two  48-position  scanivalves  and 
two  differential  pressure  tranducers  (Scanivalve  PDCR23D  with  a  zero  volume 
adapter)  were  used.  A  static-pressure  probe  was  used  to  measure  the  ambi¬ 
ent  pressure  (excluding  the  effects  of  gravity).  It  was  located  at  x  = 
-.1,  y  =  .85,  and  z  =  0.25.  This  position  was  selected  so  that  the  probe 
was  in  undisturbed  flow  and  its  wake  and  wavemaking  did  not  disturb  the 
flow  in  the  measurement  region.  All  of  the  pressure  taps  were  joined  by 
vinyl  tubing  to  eight  male  48-port  pneumatic  connectors.  Since  this 
arrangement  allows  for  the  measurement  of  376  pressures,  some  locations 
(17)  were  measured  twice  through  the  use  of  a  branch  to  insure  repeatabil¬ 
ity.  Consecutively,  two  of  the  connectors  at  a  time  were  joined  to  female 
connectors  which  were  connected  to  the  scanivalves.  Pressure  tubes  from 
the  scanivalves  were  connected  to  the  plus  side  of  the  transducers.  Tubes 
from  the  static-pressure  probe  (divided  at  a  branch)  were  connected  to  the 
negative  side  of  the  transducers.  The  tubes  from  the  branch  were  also 
connected  to  one  channel  of  each  connector  to  check  the  zero  point  of  the 
transducers  during  the  running  of  the  carriage.  The  scanivalves  were 
driven  by  solenoid  controllers  which  were  controlled  by  a  microcomputer  on 
the  carriage.  By  the  control  of  the  scanivalve  and  the  choice  of  connec¬ 
tors,  any  one  of  the  pressure  taps  on  the  hull  could  be  connected  to  the 
transducer.  Figure  5b  provides  a  block  diagram  of  the  system  used  for  the 
measurement  of  the  surface-pressure  distribution. 

Two  five-hole  pitot  probes  (modified  NPL  type),  one  for  the  port  and 
one  for  the  starboard  side  of  the  hull,  were  used  to  measure  the  direction 
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and  magnitude  of  the  velocity  in  the  hull  boundary  layer  and  wake.  The 
starboard-side  probe  is  shown  in  figure  6.  The  port-side  probe  is  similar, 
but  a  mirror  image.  Figure  7  is  a  photograph  of  the  automated  traverse 
used  to  position  the  probes.  The  probes  could  be  moved  in  three  directions 
of  a  Cartesian  coordinate  system.  The  port  and  starboard  crossplane  posi¬ 
tioning  (y-z  planes)  was  driven  by  three  stepper  motors  which  were  control¬ 
led  by  a  microcomputer  on  the  carriage.  This  enabled  multiple  measurements 
per  carriage  run.  The  axial  (x-direction)  positioning  was  achieved  by 
moving  the  entire  traverse  along  the  measurement  rails  of  the  carriage. 
The  same  static-pressure  probe  as  described  above  was  also  used.  The  leads 
from  the  pitot  probes  were  connected  by  vinyl  tubing  to  one  side  of  ten 
differential  pressure  transducers  (Kyowa  PDIOOCA).  The  tube  from  the  sta¬ 
tic-pressure  probe  was  connected  by  vinyl  tubing  to  a  branch.  Ten  pressure 
tubes  were  divided  at  the  branch  and  connected  to  the  other  side  of  the 
transducers.  A  block  diagram  of  this  arrangement  is  shown  in  figure  8a. 

For  each  experiment,  the  voltage  output  from  the  transducers  was  sam¬ 
pled,  digitized,  recorded,  and  analyzed  by  a  microcomputer  on  the  carriage. 
The  measurements  were  monitored  with  a  multi-pen  recorder  during  each  car¬ 
riage  run. 

C.  Calibracion 

The  differential  pressure  transducers  were  calibrated  using  two  water 
tanks.  One  was  moved  up  and  down  by  a  one-dimensional  traverse  while  the 
other  was  at  a  fixed  elevation,  l.e.,  the  pressure  was  measured  by  water 
head. 

Both  five-hole  pitot  probes  were  calibrated  in  the  towing  tank  using 
the  calibration  device  shown  in  figure  9.  Since  the  effect  of  Reynolds 
number  on  calibration  is  known  to  be  insignificant  for  carriage  speeds  in 
excess  of  .4m/s,  the  calibration  was  done  at  a  carriage  speed  of  Im/s.  The 
data  were  analysed  using  a  method  similar  to  Fujita  (1979)  which  is  des¬ 
cribed  below. 

The  probe-based  Cartesian  coordinate  system  (X,Y,Z)  and  other  basic 
quantities  used  for  the  data  anlysis  are  defined  in  figure  10.  Referring 
to  figure  10,  the  X-axis  is  in  the  direction  0-C  and  the  Y-  and  Z-axes  are 


in  the  T-C-B  and  P-C-S  planes,  respectively.  The  origin  is  at  the  center 
of  curvature  of  the  tip  face.  The  probe  holes  are  designated  C,  T,  B,  P, 
and  S  and  the  corresponding  sensed  pressures  are  Hg,  Hp,  and  Hg , 
respectively.  The  calibration  velocity  vector  is  V  and  its  projections  in 
the  X-Y  and  Y-Z  planes  are  and  V^,  respectively.  The  angles  between  Vj^ 
and  Vy  and  the  X-axis  are  6^  and  g^,  respectively.  The  angle  between  V  and 
Vy  is  9j^.  As  described  above,  etc.,  were  measured  as  the  differential 
pressures  relative  to  the  static-pressure  tube.  The  probes  were  calibrated 
in  the  range  -45°  9^^  ^  ^5°  and  -45°  8^  ^  ^5°  with  five-degree  steps 
in  9j^  and  8^.  The  following  calibration  coefficients  were  then  calculated: 


M(9j^,6^)  =  (4H^-  H^-  Hg-  Hp-  H^j/CV^/Zg)  (II-l) 
K(9^,6^)  =  (H^-  H3)/(4H^-  H^-  Hg-  Hp-  Hg)  (II-2) 
L(9h.Bv)  =  (Hg-  Hp)/(4H^-  H^-  Hg-  Hp-  Hg)  (II-3) 
P(0j^,6^)  =  H^/(v2/2g)  (II-4) 


2 

where  g  is  the  gravitational  acceleration,  V  =  V  V,  and  tan  9^^  = 
cos  g^tangj^.  These  coefficients  are  shown  in  figure  11. 

The  procedure  for  evaluating  the  mean-velocity  components  (u,v,w)  and 
pressure  p  from  the  local  probe  measurements  and  the  above  calibration 
coefficients  is  as  follows.  The  probe  is  positioned  such  that  local  probe 
coordinates  (X,Y,Z)  are  parallel  to  the  global  hull  coordinates  (x,y,z). 
Designating  local  measured  values  with  a  ’,  the  following  quantities  are 
determined: 

M’  =  4H^  -  H^  -  H^  -  H^  -  H^  (II-5) 
K'  =  (H^  -  H^)  /  (4H^  -  H^  -  H^  -  H^  -  H^)  (11-6) 
L’  =  (H^  -  Hp  /  (4H^  -  H^  -  H^  -  H^  -  H^  (11-7) 
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0'  and  6'  are  determined  from  K’  and  L'  using  the  calibration  coefficients 
n  V 

shown  in  figures  11a  and  11b.  These  are  then  used  to  determine  M  and  P, 
respectively,  from  figures  11c  and  He,  or  lid  and  Ilf.  As  a  result,  the 
following  are  obtained: 


8,  '  =  tan  ^  (tan0 '/cos6 ' ) 
h  h  V 


(II-8) 


V’  -  /2gM'/M 


(II-9) 


V*  =  V'cos8’//l“Sin^B’sln^B,' 

h  V  V  h 


(II-IO) 


V  =  V'cosBu/*^l~8in^6'sin^B,' 

V  h  V  h 


C  =  2(p  -  p^)/pU^  =  (HI  -  PV^/2g)/(U^/2g) 

P  O  w 


(11-11) 


(111-12) 


where  p^  and  U  are  the  pressure  and  velocity  in  the  undisturbed  uniform 
stream.  Finally,  the  velocity  components  are  determined  as: 


Vi;cosB^ 


(11-13) 


V'sins; 


(11-14) 


w  =  V'sinB' 

V  V 


(11-15) 


Figure  8b  provides  a  block  diagram  of  the  overall  data  analysis  pro¬ 
cedure.  A  computer  program  was  used  for  data  analysis  in  which  values  at 
desired  points  are  obtained  by  Lagrange  Interpolation. 


D.  Data-Acquisitlon  System 

The  data-acquisition  system  is  a  microcomputer  PC9801  VM2  with  two  8 
channel  A-D  converter  boards  and  a  digital  10  board  which  controlled  the 
stepper  motors  and  scanivalves.  The  pressures  sensed  from  either  the 
probes  or  pressure  taps  were  converted  to  voltage  by  the  transducers  and 
then  filtered  by  a  low-pass  filter  and  sampled  through  the  A-D  converter. 
A  sampling  frequency  of  SHz  was  used  in  all  the  experiments. 
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E.  ExperlaeaCal  Procedures 

Six  types  of  measurements  were  made;  resistance,  propeller  open-water, 
self-propulsion,  wave-profiles,  surface-pressure  distribution,  and  mean- 
velocity  and  pressure  field.  In  all  cases,  a  time  interval  of  about  15 
minutes  between  carriage  runs  (i.e.,  maximum  of  four  runs  per  hour)  was 
necessary  in  order  for  the  fluid  motion  induced  by  the  previous  run  to  be 
sufficiently  damped.  All  measurements  are  for  the  full-load  condition. 
First,  resistance  and  self-propulsion  tests  as  well  as  wave-profile  mea¬ 
surements  were  made.  Based  on  these,  the  conditions  for  the  subsequent 
detailed  measurements  were  selected. 

The  resistance  and  propulsion  tests  were  carried  out  following  stan¬ 
dard  towing-tank  procedures.  Both  tests  were  for  the  model-free  condition 
(i.e.,  model  was  free  to  sink  and  trim).  The  dynamometers  were  calibrated 
before  and  after  the  measurements  in  all  experiments.  Force  measurements 
were  made  for  about  20  seconds  after  the  carriage  attained  steady  speed. 
An  automated  clamp  was  used  to  secure  the  hull  while  the  carriage  was 
accelerating  or  decelerating.  Measurements  were  performed  for  the  following 
conditions:  carriage  speed,  *5  <_ U  ^  1.7m/s;  Froude  number,  .08  ^  Fr  ^  *27; 
and  Reynolds  number,  1.6  x  10^  _<_  Re  _<_  5.4  x  10^. 

The  propeller  open-water  test  was  performed  with  the  number  of  propel- 

2  S 

ler  revolutions  constant,  i.e.,  n  =  lOrps  (nDp/v  *  1.7  x  10^).  The  car¬ 
riage  speed  range  was  0  ^  1.55m/s  (0  ^  ^  1.1).  Self-propulsion  tests 
were  performed  for  six  speeds  U  *  (.7,  .8,  .9,  1,  1.1,  1.2m/s).  The  corre¬ 
sponding  Froude  numbers  are  Fr  =  (.112,  .128,  .144,  .16,  .176,  .192).  For 
each  speed,  thrust,  torque,  and  resistance  were  measured  with  the  load 
varying  from  zero  to  the  model  self-propulsion  point. 

The  wave  profiles  were  recorded  only  for  the  model-fixed  condition 
(i.e.,  model  was  fixed  at  the  design  draft).  This  was  done  photograph¬ 
ically  using  both  35mm  and  video  cameras,  and  normal  processing. 

For  the  detailed  measurements,  the  speed  was  selected  as  U  =  Im/s  in 
order  to  minimize  free-surface  effects  and  to  aid  in  maintaining  measure¬ 
ment  accuracy.  The  number  of  propeller  revolutions  for  the  with-propeller 
condition  was  7,8rps,  and  the  corresponding  thrust  and  torque  were  .0667!! 


and  .0017Nm,  respectively.  These  are  all  close  to  the  model  self-propul¬ 
sion  point  for  this  model  speed.  Under  these  conditions,  the  propeller 
influence  is  large.  All  the  measurements  were  made  for  the  model-fixed 
condition.  The  important  nondlmensional  parameters  have  the  following 
values  (see  List  of  Symbols  for  definitions): 

Fr  -  U//gL  =  .16 

Re  -  UL/v  ~  3.94  x  10^ 

C,j  -  2T/pU^Tr(Rp-  rj)  =■  .803 

(11-16) 

=  T/pn^Dp^  =  .234 
Kq  »  Q/pn^Dp^  =  .0411 
J  »  U^/nDp  =  .654 

For  the  surface-pressure  distribution,  three  measurements  were  made 
per  carriage  run.  The  first  measurement  was  used  to  check  the  zero  point 
and  the  following  two  to  measure  pressures  at  four  points  by  using  the  two 

transducers.  The  change  was  carried  out  every  20m  (i.e.,  20s)  with  a  pho¬ 

to-switch  system.  Measurements  were  made  at  nine  axial  stations  and  along 
three  waterlines  for  both  the  with-  and  without-propeller  conditions,  on 
both  the  port  and  starboard  sides  of  the  hull. 

For  the  mean-velocity  and  pressure  fields,  four  measurements  were  car¬ 
ried  out  per  carriage  run  (i.e.  ,  eight  points  were  measured  by  use  of  the 
port  and  starboard  probes).  For  the  Im/s  condition,  the  carriage  ran 
steadily  for  about  65  seconds.  It  took  about  one  second  to  traverse  the 
probe  from  one  point  to  the  next,  and  the  probe  response  time  was  about  two 
seconds.  Therefore,  ten  seconds  were  used  to  make  measurements  at  each 
point,  with  the  last  six  seconds  used  to  obtain  averages.  The  traverse  was 

carried  out  every  15m  by  a  photo-switch  on  the  carriage  and  shutter  plates 

on  the  cat  walk.  The  probes  usually  had  small  setting  angles  (1  to  2°), 
and  in  some  cases  fairly  large  preset  angles  (around  10°)  so  that  measure- 
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merits  could  be  performed  close  to  the  hull.  The  flow  angles  were  corrected 
based  on  measurements  in  uniform  flow.  The  mean-velocity  and  pressu’'e 
field  measurements  were  also  performed  '“»•  both  conditions,  i.e.,  with  and 
without  propeller,  and  on  both  sides  of  the  hull,  for  up  to  thirteen  axial 
stations  (.5  _<  x  _<  1.1).  For  stations  far  upstream  of  the  propeller  and, 
of  course,  the  propeller  plane  itself,  the  measurements  were  made  for  the 
without-propeller  condition  only.  Measurements  at  about  400  -  700  points 
were  made  at  each  transverse  (x  =  constant)  section.  The  locations  for 
both  the  surface-pressure  distribution  and  mean-velocity  and  pressure  field 
measurements  are  summarized  in  Table  2. 

F.  Experimental  Oncertalnty 

The  accuracy  of  the  mean-velocity  measurements  is  estimated  to  be 
within  1.5%  for  the  magnitude  and  1  degree  for  the  direction.  The  accuracy 
of  the  measurement  of  the  pressure  coefficient  Cp  is  estimated  to  be  within 
±  .01  and  ±  .05  for  the  surface  and  field  values,  respectively. 

III.  RESULTS 

In  the  following,  the  detailed  experimental  results  are  presented  and 
discussed  to  point  out  the  essential  features  of  the  flow  for  both  the 
with-  and  without-propeller  conditions.  In  most  cases,  the  results  for  the 
without-propeller  condition  are  discussed  first,  followed  by  those  for  the 
with-propeller  condition.  Although  not  discussed  in  the  Introduction,  in 
comparison  to  the  situation  for  propeller-hull  interaction,  a  considerable 
amount  of  detailed  experimental  Information  is  available  for  the  mean-flow 
over  the  stern  and  in  the  near  wake  for  bare  ship  hulls  either  for  double 
bodies  or  at  low  Fr.  Moat  of  the  experiments  are  for  merchant  ships,  fair¬ 
ly  similar  to  the  Series  60  Cg  =  .6,  but  with  larger  block  coefficients 
(Cg  ~  .8).  Several  interpretations  have  been  given  to  the  vortical  flow 
patterns  (and  attendant  complicated  velocity  contours)  associated  with  the 
afterbody  flow  (i.e.,  the  flow  over  the  stern  and  in  the  near  wake)  for 
these  hull  forms.  This  is  exemplified  by  the  range  of  terminology  used  to 
describe  them,  e.g.,  bilge  vortices,  three-dimensional  separation,  and 
longitudinal  vorticity.  For  a  recent  review  of  this  topic,  see  Patel 
(1988).  Below,  we  shall  simply  point  out  these  features  with  the  primary 
emphasis  on  their  modification  due  to  the  action  of  the  propeller. 
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A.  Resistance  and  Propulsion  Tests 

Resistance  tests  were  performed  to  provide  data  for  evaluating  the 
thrust-deduction  factor.  Also,  as  will  now  be  discussed,  it  was  of  inter¬ 
est  to  compare  the  present  results  with  those  of  the  CEP  (see  Chapter  I). 

Figure  12  shows  a  comparison  of  the  present  results  for  the  total- 
resistance  coefficient  C.^  with  one  of  those  obtained  in  the  CEP,  namely, 
Akishima  Laboratories,  MITSUI  ZOSEN  (ALM).  Although  the  trends  are  very 
similar,  the  present  values  are  slightly  smaller.  However,  of  the  four 
sets  of  results  of  the  CEP  for  4m  models,  the  ALM  results  indicated  the 
largest  values. 

The  residuary-resistance  coefficient  Cp^  is  shown  in  figure  13  along 
with  the  three  mean  Cp  curves  for  the  different  model-size  groups  of  the 
CEP.  Although  the  present  results  are  in  agreement  with  those  for  the 
larger  model-size  group,  they  show  somewhat  larger  humps  and  hollows  than 
the  mean  curves.  This  is  due  to  the  averaging  technique  used  to  define  the 
mean  curves  since  similar  humps  and  hollows  at  nearly  the  same  Fr  are  evi¬ 
dent  in  the  individual  curves,  e.g. ,  the  10m  (NK.K)  and  2,5m  (UT)  model 
results  in  the  original  reference  (ITTC,  1987), 

Based  on  the  above  comparisons,  it  is  concluded  that  the  present  re¬ 
sults  are  in  close  agreement  with  those  obtained  in  the  CEP,  in  spite  of 
the  small  modification  made  to  the  stern  geometry  to  attach  the  propeller 
(see  Section  11. A  and  figure  3). 

The  propeller  open-water  test  was  performed  to  investigate  the  propel¬ 
ler  performance  in  uniform  flow  and  to  provide  data  for  evaluating  the 
effective-wake  factor  and  relative  rotative  efficiency  based  on  the  thrust- 
identity  method.  The  results  are  shown  by  open  symbols  in  standard  format 
in  figure  14.  As  expected,  the  value  of  J  at  which  ICj,  =  0  is  a  little 
larger  than  the  geometrical  pitch  ratio.  Also,  note  the  rather  wide  range 
of  J  for  which  the  propeller  efficiency  is  large.  The  present  results  are 
in  close  agreement  (figure  not  shown)  with  those  of  the  MAU  methodical 
series  tests  of  Tsuchida  et  al.  (1958). 
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The  self-propulsion  test  results  are  shown  in  figure  15.  Figure  15a 
shows  the  results  for  all  six  speeds  tested,  whereas  figure  15b  is  for  U  * 
In/s  which  is  the  speed  for  the  detailed  measurements.  As  expected  from 
physical  considerations,  the  towing  force  R  decreases  with  increasing 
thrust  T.  The  decrease  of  R  is  smaller  than  the  increase  of  T.  The  R  vs. 
T  relationship  is  almost  linear  with  the  same  slope  (about  -.9)  for  each 
U.  Therefore,  if  this  slope  is  used  to  define  the  thrust-deduction  factor, 
it  would  be  nearly  constant  for  all  load  conditions  and  Fr.  However,  the 
usual  definition,  i.e. 

(1-t)  =  (Rq  -  R)/T  (III-l) 

where  R^  is  the  without-propeller  resistance,  is  used  below.  Figure  15b 
also  shows  Rq  and  (T  +  R)  vs.  T.  It  is  seen  that  the  interaction  between 
the  propeller  and  hull  is  weak  for  T  =  0  (i.e.,  R(T  =  0)  is  slightly  larger 
than  Rq)  and  increases  gradually  with  T. 

The  usual  self-propulsion  factors  are  shown  vs.  thrust-loading  coeffi¬ 
cient  Cj  for  each  U  in  figure  16.  As  already  mentioned,  the  effective-wake 
faCor  (l-w)  and  relative  rotative  efficiency  (nj^)  were  determined  based  on 
the  thrust-identity  method  (e.g.,  Todd,  1967),  i.e. 


(l-w)  - 

''R  ■ 


(III-2) 


where  is  the  speed  of  advance  and  Qq  the  torque  in  open  water,  i.e.,  the 
speed  and  torque  in  uniform  flow  at  which  the  propeller  would  produce  the 
same  thrust  at  the  same  revolutions  per  second  as  that  measured  behind  the 
model.  The  effective-wake  factor  (l-w),  shown  in  figure  16a,  increases 
gradually  with  C.^  due  to  the  propeller-induced  flow  contraction  and  axial 
acceleration.  This  is  consistent  with  previous  work  in  which  it  is  also 
shown  that  the  increase  can  be  significant  for  hulls  with  large  block  coef¬ 
ficients  (e.g.,  Nagamatsu  and  Sasajima,  1975).  Note  that  the  nominal 
volume  mean  velocity  for  U  =  Im/s  is  .66.  The  thrust-deduction  factor  (1  - 
t),  shown  in  figure  16b,  also  increases  gradually  with  Cq<,  in  this  case. 
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due  to  Che  resistance  increase  for  T  *  0.  Lastly,  the  relative-rotative 
efficiency  n|^,  shown  in  figure  16c,  is  nearly  constant,  but  shows  some 
scatter  for  small  due  to  difficulties  in  measuring  small  values  of 

torque.  For  the  low  Fr  range  of  the  present  self -propulsion  tests  (.1  <  Fr 
<  .2),  the  self-propulsion  factors  have  almost  constant  values  if  plotted 
vs.  Fr  for  fixed  C^,  which  is  the  usual  format,  e.g.,  at  the  model  self¬ 
propulsion  point  C.p  =  .8,  (1  -  w)  “  “  .74-. 75,  (1  -  t)  *  .86-. 87,  and 

=  1-1.02  for  all  Fr.  However,  based  on  previous  work  (e.g.,  Yamazaki 
and  Nakatake,  1984),  these  factors  are  expected  to  vary  with  Fr  for  larger 
Fr. 


B.  Wave  Profiles 

The  wave  profiles  at  the  hull  were  measured  using  both  35mm  and  video 
cameras.  The  profiles  for  Fr  =  (.3,  .25,  .16)  are  shown  in  figure  17, 
including  for  the  former  two  cases,  comparisons  with  results  from  the  CEP. 
Although  there  are  some  differences  in  the  experimental  conditions,  the 
present  results  are  in  close  agreement  with  those  of  the  CEP.  For  Fr  “ 
.16,  which  is  the  condition  for  the  detailed  measurements,  results  are 
shown  for  both  the  with-  and  without-propeller  conditions.  It  is  seen  that 
the  propeller  influence  is  negligible  for  x  ^  .9.  For  .9  <  x  <  .98,  the 
wave  elevation  with  propeller  is  less  than  without,  but  for  .98  £  x _<  1  the 
reverse  holds  true.  As  will  be  shown  below,  this  correlates  with  the  pres¬ 
sure  measurements  and  is  associated  with  the  propeller  thrust. 

C.  Surface-Pressure  Distributions 

The  surface-pressure  distributions  are  shown  in  figures  18  through 
20.  The  pressure-tap  locations  are  shown  in  figure  5a.  Figure  18  shows 
the  pressure  contours  for  both  the  with-  and  without-propeller  conditions 
as  well  as  for  the  difference  between  the  two  conditions.  The  results 
presented  are  for  the  port  side  only  since,  as  will  be  shown  below,  the 
distributions  are  nearly  symmetric.  For  the  without-propeller  condition, 
it  is  seen  that  the  pressure  is  relatively  low  near  the  keel  and  the  pres¬ 
sure  recovery  is  larger  for  waterlines  near  the  free  surface.  For  the 
with-propeller  condition,  the  expected  pressure  decrease  is  observed  near 
the  propeller,  except  above  the  propeller,  where  there  is  an  increase.  The 
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pressure  decrease  is  larger  in  the  lower  part  of  the  propeller  disk  than 
the  upper,  where  the  boundary  layer  is  thick.  At  large  distances  upstream 
from  the  propeller  (x  <  .95),  the  pressure-difference  contours  are  uniform 
with  depth,  whereas  for  smaller  distances,  they  are  curved  towards  the 
outer  propcxler  radii.  The  influence  of  the  propeller  is  restricted  to 
X  >  .9  (which  is  about  two  diameters  upstream  of  the  propeller). 

The  girthwise  and  streamwise  variations  (figures  19  and  20,  respec¬ 
tively)  show  the  same  trends  as  just  described.  Figure  19  shows  that  the 
magnitude  of  the  pressure  decrease  is  almost  the  same  from  the  keel  to  the 
free  surface  for  x  <  .95,  but  for  .95  ^  1,  there  is  a  large  girthwise 
variation  indicating  a  pressure  decrease  except  above  the  propeller.  Fig¬ 
ure  20  shows  that  for  the  without-propeller  condition,  the  port  and  star¬ 
board  pressure  variations  are  very  nearly  symmetric.  The  asymmetry  is,  no 
doubt,  due  to  slight  variations  in  the  geometry  of  the  pressure  taps.  The 
with-propeller  condition  results  are  less  symmetric,  but  the  differences 
between  the  port  and  starboard  trends  are  difficult  to  discern.  Near  the 
propeller  plane,  the  magnitude  of  the  pressure  decrease  clearly  depends  on 
the  propeller  spanwlse  loading  distribution  (i.e.,  (z-Zjj)/Rp).  Also,  the  z 
=  .008  results  (figure  20a)  are  quite  similar  to  the  wave  elevation  (cf. 
figure  17). 

D.  Mean-Velocity  and  Pressure  Fields 

The  mean-velocity  and  pressure  fields  for  the  with-  and  without- 
propeller  conditions  are  shown  in  figures  21  through  23.  Figure  21  shows 
the  mean-velocity  vectors  projected  onto  the  crossplane  and  the  contours  of 
axial  velocity.  Figures  22  and  23  show  the  axial  vorticity,  =  (v^  - 
Wy),  and  pressure  contours,  respectively. 

Considering  the  without-propeller  condition  first,  it  is  seen  from 
figure  21a  that,  at  the  first  measurement  station,  x  *  .5,  the  flow 
diverges  off  the  bilge  such  that  the  boundary  layer  is  thicker  near  the 
centerplane  and  the  free  surface  than  in  the  bilge  region.  However,  at  x  = 
.6  (figure  21b),  which  may  be  regarded  as  the  beginning  of  the  afterbody 
flow,  the  flow  diverges  off  the  centerplane  such  that  the  boundary-layer 
thickness  is  nearly  uniform  around  the  entire  girth.  For  both  these  sta- 
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tions  (x  “  .5,  .6),  the  crossplane  flow  Is  very  small.  It  appears  that, 
for  this  particular  experiment,  x  =  .6  may  be  an  ideal  station  to  initiate 
stern-flow  calculations  with  initial  conditions  based  on  simple  two- 
dimensional  boundary-layer  correlations.  Subsequently,  at  x  *  .7  and  .8 
(figures  21c, d),  the  flow  begins  to  exhibit  features  that  are  well  known 
for  this  type  of  hull  form,  features  which  become  quite  evident  in  the 
stern  region,  x  =  .9,  .95,  .975,  and  .98125  (figures  21e-h).  That  is,  the 
diminishing  cross-section  of  the  hull  at  the  stern  and  the  associated  con¬ 
vergence  of  the  inviscld  streamlines  leads  to  a  thickening  of  the  region  of 
viscous  flow,  except  near  the  keel  where  there  is  a  thinning  due  to  flow 
divergence  off  of  the  centerplane.  Generally,  the  crossplane  flow  is 
directed  upwards  and  towards  the  hull  centerplane.  Also,  as  will  be  dis¬ 
cussed  further  below,  quite  apparent  in  the  crossplane  velocity  vectors  is 
the  presence  of  a  weak  longitudinal  vortex.  The  pronounced  bulge  in  bound- 
ary-layer  thickness  near  the  region  of  maximum  hull  concavity  and  the 
rather  complicated  flow  in  the  vicinity  of  the  stern  Lube  are  particularly 
noteworthy.  Finally,  in  the  region  .9875  x  _<  1,1  (figures  211-m),  the 
near  wake  shows  relatively  fast  decay  of  the  transverse  velocity  components 
and  the  initial  stages  of  wake  recovery.  There  is  a  gradual  increase  in 
the  wake  centerplane  velocity  and  diffusion  of  the  wake.  Mote  the  change 
in  the  shape  of  the  wake  as  it  evolves,  indicating  a  relatively  slower 
recovery  near  the  free  surface  and  wake  centerplane.  Some  of  these  trends 
are  due  to  the  downward  convection  associated  with  the  weak  longitudinal 
vortex  and  the  general  upward  motion  of  the  external  inviscid  flow  in  this 
region. 

For  the  with-propeller  condition,  at  x  =  .9  (figure  21e),  the  influ¬ 
ence  of  the  propeller  is  negligible,  although  there  is  a  slight  decrease  in 
boundary-layer  thickness.  Subsequently,  at  x  =  .95,  .975,  and  .98125  (fig¬ 
ures  21f-h),  the  propeller  influence  becomes  increasingly  apparent  as  the 
propeller  plane  (x  =  .9875)  is  approached,  especially  in  the  axial  velocity 
contours  which  exhibit  increased  velocity  near  the  hull  due  to  the  combined 
effects  of  the  propeller-induced  flow  contraction  and  axial  acceleration. 
Also,  evident  in  the  crossplane  vectors  is  a  slight  turning  of  the  flow, 
i.e.  ,  the  crossplane  vectors  are  turned  towards  the  propeller  axis  and 
increase  slightly  in  magnitude  due  to  a  small,  but  discernible,  propeller- 
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induced  radial  velocity.  The  vortical  flow  is  still  present,  but  appears 
to  be  compressed  towards  the  hull.  Note  that  the  flow  upstream  of  the 
propeller  is  nearly  symmetric  with  respect  to  the  hull  centerplane.  The 
flow  in  the  propeller  slipstream,  x  =  1,  1.01875,  1.05,  and  1.1  (figures 
21j-m),  is,  of  course,  completely  altered  due  to  the  action  of  the  propel¬ 
ler,  indicating  characteristics  which  are  similar  to  a  swirling  jet  and  a 
complex  interaction  between  the  hull  boundary  layer  and  wake  and  the  pro¬ 
peller-induced  flow.  Quite  obvious  is  the  propeller-induced  flow  asymmetry 
within  the  propeller  slipstream  with  respect  to  the  hull  centerplane, 
including  a  very  significant  drift  of  the  wake  centerplane  off  the  hull 
centerplane  towards  the  port  side  of  the  hull.  Outside  the  propeller  slip¬ 
stream,  the  flow  is  nearly  symmetric  with  respect  to  the  hull  center- 
plane.  The  swirl  velocity  is  maximum  just  downstream  of  the  propeller 

(figure  21j)  and  then  decays  relatively  rapidly  in  the  near  wake  (figures 
21k-m).  The  swirl  profiles  reflect  the  complex  interaction  between  the 
propeller  and  the  hull  boundary  layer  and  wake,  and  cannot  be  explained 
simply  as  a  superposition  of  the  bare-hull  and  propeller-induced  flows. 
Figures  21k-m  also  suggest  the  occurrence  of  a  secondary  vortex  on  the 
starboard  side  of  the  hull  near  the  free  surface.  The  axial  velocity  con¬ 
tours  indicate  that  the  slipstream  flow  initially  accelerates  (figures  21j- 
l)  and  then  undergoes  a  rapid  deceleration  and  gradual  diffusion  (figure 
21m).  Considerably  higher  velocities  are  found  on  the  starboard  than  on 
the  port  side  of  the  hull.  The  concentration  of  the  axial  velocity  con¬ 
tours  in  figures  21j-Ji,  are,  no  doubt,  the  trace  of  the  blade-tip  vortices 
and  its  associated  vortex  sheet  in  the  near  wake  (figures  21j,k)  which  then 
diffuses  with  downstream  distance  (figures  212,, m). 

Next,  we  consider  the  characteristics  of  the  axial  vorticlty  contours 
shown  in  figure  22.  These  were  determined  by  differentiation  of  the  cross¬ 
plane  velocity  field,  made  possible  by  the  closely-spaced  Cartesian  grid 
used  in  the  measurements.  For  the  without-propeller  condition,  the  axial 
vorticlty  is  small  at  the  upstream  stations  (.5  x  _<  .8)  where  the  cross¬ 
plane  flow  is  also  small,  and  therefore,  the  flow  is  nearly  two-dimensional 
(figures  22a-d).  Subsequently  (x  =  .9  and  .95),  there  is  a  significant  in¬ 
crease  in  the  axial  vorticlty.  This  is  primarily  due  to  the  developing 
three-dimensionality  of  the  flow  and  the  associated  rotation  of  the  bound- 
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ary-layer  vortlcity  towards  the  axial  direction  (figures  22e,f).  Finally, 
in  the  stern  and  near  wake  regions  (.975  ^  x  ^  1.1),  the  axial  vorticity 
exhibits  the  expected  pattern  associated  with  the  rapid  boundary-layer 
thickening  and  weak  longitudinal  vortex  typical  for  this  type  of  hull  form 
(figures  22g-m).  Note  the  rather  slow  decay  rate  of  the  axial  vorticity 
which  can  be  correlated  with,  but  is  much  slower  than,  that  of  the  cross¬ 
plane  flow.  For  the  with-propeller  condition,  at  x  =  .9  and  .95  (figures 
22e,f),  the  influence  of  the  propeller  is  negligible.  At  x  =  .975  and 
.98125  (figures  22g,h),  the  axial  vorticity  shows  an  increase  and  a  change 
in  shape  near  the  hull  (i.e.,  the  contours  are  compressed  towards  the  hull 
centerplane)  due  to  the  propeller-induced  flow  contraction  and  axial  accel¬ 
eration.  In  the  propeller  slipstream,  x  =  1,  1.01875,  1.05,  and  1.1  (fig¬ 
ures  22j-m),  the  axial  vorticity  is  very  large.  The  regions  of  large  nega¬ 
tive  and  positive  vorticity  correspond,  respectively,  to  the  vortices  from 
the  blade-tips  and  the  propeller  hub.  The  flow  asymmetry  within  the  pro¬ 
peller  disk  with  respect  to  the  hull  centerplane  is  quite  apparent.  The 
decay  of  axial  vorticity  with  downstream  distance  again  correlates  with  the 
crossplane  velocity  field  described  earlier,  but  with  a  considerably  slower 
rate.  Also,  note  the  presence  of  both  the  primary  and  secondary  propeller- 
induced  vortices. 

Lastly  for  this  chapter,  we  consider  the  characteristics  of  the  pres¬ 
sure  field  shown  in  figure  23.  From  the  outset  we  should  point  out  that 
the  measurement  of  pressure  in  the  flow  field  is  difficult  and  the  accuracy 
is  limited  (see  Section  II. F).  As  a  result,  the  port  and  starboard 
results,  even  for  the  wlchout-propeller  condition,  exhibit  some  anomalous 
differences  and  scatter.  In  general,  the  trends  for  both  the  with-  and 
without-propeller  conditions  are  consistent  with  the  surface-pressure  mea¬ 
surements  (figures  18  through  20),  but  with  some  differences  in  the  pres¬ 
sure  magnitude,  especially  for  the  port  side  of  the  hull,  due  to  the  dif¬ 
ferences  in  the  measurement  technique  and  analysis  methods.  For  the  with¬ 
out-propeller  condition,  the  pressure  is  largest  near  the  hull  and  grad¬ 
ually  decreases  across  the  viscous-  and  Inviscld-f low  regions.  The  cross¬ 
plane  variations  are  surprisingly  uniform  in  view  of  the  crossplane  veloc¬ 
ity  field,  although  relatively  lower  pressures  are  found  in  the  core  region 
of  the  vortical  flow.  For  the  with-propeller  condition,  at  x  *  .9  (figure 
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23e) ,  the  influence  of  the  propeller  Is  negliglible.  Subsequently,  at  x  = 
.95,  .975,  and  .98125  (figures  23f-h),  the  influence  becomes  increasingly 
apparent  as  the  propeller  plane  is  approached.  There  is  a  large  drop  in 
pressure  associated  with  the  propeller  thrust  and  the  shape  of  the  pressure 
contours  is  changed  remarkably  due  to  the  propeller-induced  flow  contrac¬ 
tion  and  axial  acceleration,  the  pressure  contours  being  nearly  circular 
within  the  propeller  disk.  The  pressure  contours  upstream  of  the  propeller 
are  nearly  symmetric  with  respect  to  the  hull  centerplane.  Immediately 
downstream  of  the  propeller,  at  x  =  1  (figure  23j),  there  is  a  large 
increase  in  pressure,  associated  with  the  propeller  thrust,  except  at  the 
center  of  the  slipstream  where  there  is  a  low  pressure  region  associated 
with  the  hub  vortex.  At  this  station,  the  pressure  contours  are  asymmetric 
with  respect  to  the  hull  centerplane.  Considerably  larger  pressures  are 
found  on  the  starboard  than  on  the  port  side  of  the  hull.  Finally,  at  x  = 
1.05,  1.01875,  and  1.1  (figures  23k-m),  there  is  a  rapid  decrease  in  pres¬ 
sure.  This  is  consistent  with  the  previously  mentioned  initial  axial 
velocity  increase  in  this  region.  The  pressure  field  changes  rapidly  such 
that  the  contours  become  circular  and  nearly  symmetric  with  respect  to  the 
wake  centerplane.  The  continued  presence  of  the  low  pressure  region  near 
the  center  of  the  slipstream  indicates  the  hub  vortex.  It  is  interesting 
to  observe  that,  unlike  the  axial  velocity  contours,  the  pressure  contours 
have  a  similar  shape  (i.e.,  circular)  within  the  propeller  disk  both 
upstream  and  downstream  (except  immediately  downstream)  of  the  propeller. 

IV.  PROPELLER-HOLL  INTERACTION 

The  foregoing  discussion  of  the  experimental  results  clearly  indicates 
the  complexity  of  the  interaction  between  the  propeller-induced  flow  and 
the  hull  boundary  layer  and  wake.  In  this  section,  the  results  are  anal¬ 
yzed  further  to  assess  the  nature  of  this  interaction.  To  this  end,  use  is 
made  of  the  Kerwln  and  Lee  (1978)  propeller-performance  program  with  both 
nominal  and  effective  Inflows.  The  effective  inflow  was  calculated  using 
the  method  of  Toda  et  al.  (1987),  who  extended  the  method  of  Huang  and 
associates  (1976,  1980)  for  axisymmetric  bodies  to  three-dimensional 
bodies.  Since  the  method  of  Toda  et  al.  (1987)  is  of  general  interest,  an 
English  translation  has  been  provided  as  Appendix  I. 
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Some  modifications  of  the  propeller-performance  program  were  required 
in  order  to  perform  the  present  calculations.  These  included  changes  to 
handle  the  MAU  propeller  geometry  and,  as  discussed  further  below,  to  eval¬ 
uate  the  time-averaged  angular  and  radial  variation  of  thrust  and  torque 
for  three-dimensional  nonuniforn  inflow.  In  the  discussions  to  follow,  a 
cylindrical  coordinate  system  (x,r,9)  is  adopted  in  which  x  coincides  with 
the  propeller  shaft  and  is  positive  downstream,  r  is  the  radial  distance 
from  the  x-axis,  and  9  is  the  angle  measured  from  top  dead  center  and 
positive  clockwise  looking  upstream  (figure  4b).  The  mean-velocity  compon¬ 
ents  in  the  directions  of  the  coordinate  axes  are  denoted  by  (V^,Vj.,Vg). 

Before  discussing  propeller-hull  interaction,  it  is  helpful  to  first 
examine  the  propeller  performance.  Initial  calculations  were  made  for 
uniform  Inflow,  i.e.  ,  the  open-water  condition.  Referring  to  figure  14,  it 
is  seen  that  the  calculations  are  in  close  agreement  with  the  present 
experimental  results.  Also,  the  calculated  spanwise  circulation  distribu¬ 
tion  r  shown  in  figure  24b  is  typical,  including  the  maximum  value  at 

r /Rp  ”  . 7 . 

Next,  consideration  is  given  to  the  results  for  the  nominal  inflow. 
The  nominal  inflow  was  obtained  from  the  measurements  in  the  propeller 
plane  (x  =  .9875)  for  the  without-propeller  condition  (figure  21i).  First, 
the  velocity  components  were  transformed  from  the  Cartesian  (x,y,z)  to  the 
cylindrical  (x,r,9)  coordinates,  and  then  Interpolated  onto  a  fixed  set  of 
radii  and  expanded  in  a  Fourier  series  in  9.  The  axial  and  tangential 
V  velocity  components  obtained  for  r/Rp  =  .7  are  shown  in  figure  25.  The 
circumferential-average  axial  velocity  is  shown  vs.  r/Rp  in  figure  24a. 
Figure  24  shows  both  the  circumferential-average  (figure  24b)  and  angular 
variation  (figure  24c)  of  the  spanwise  circulation  distribution  F.  It  is 
clear  from  figure  24b  that,  for  nonuniform  inflow  (i.e.,  tne  nominal 
inflow),  the  circumferential-average  F  changes  significantly  compared  to 
the  uniform-inflow  distribution,  F  being  considerably  larger  near  the  root 
and  slightly  smaller  near  the  tip  than  that  for  uniform  inflow.  With  non- 
uniform  inflow,  F^^^  is  larger  and  occurs  at  r/Rp  =  .6.  These  changes  are 
due  to  the  Increase  near  the  root  and  decrease  near  the  tip  of  the  circum¬ 
ferential-average  section  angle  of  attack  for  the  nonuniform-inflow  condi¬ 
tion.  Note  that  J  =  .581  based  on  the  nominal  volume  mean  velocity.  The 
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angular  variation  of  T  shown  in  figure  24c  can  also  be  correlated  with 
variations  in  the  section  angle  of  attack: 


a  = 


(IV-1) 


with 

V  +  V 

6,  =  tan-1  ^  (IV-2) 

I  2rrn  -  Vp^ 

In  (IV-1),  <ji  is  the  geometric-pitch  angle  and  the  hydrodynamic-pitch 
angle  defined  by  (IV-2),  in  which  Vpx  and  Vpg  are  the  axial  and  tangential 
components,  respectively,  of  the  propeller-induced  velocity,  which  depend 

upon  the  inflow  velocity  (V  and  V  ).  Thus,  we  see  that  r  is  large  for 
0°  <  9p  <  90®  and  9p  =  -  150®  mainly  due  to  the  effects  of  the  nominal 
axial  velocity  V^^,  but  also  due  to  the  influence  of  the  nominal  tangential 
velocity  V  (i.e.,  influence  of  the  crossplane  flow  which  is  directed 

upwards  and  towards  the  centerplane)  for  9p  s'  ±  90°. 

For  a  three-dimensional  nonuniform  inflow,  the  total  unsteady  blade 
thrust  K.J  includes  contributions  from  both  circulation  and  added  mass 

(Kjg).  These  are  shown  vs.  9p  in  figure  25.  Results  are  shown  both  with 
and  without  the  nominal  tangential  velocity  V  .  Also  shown  is  the  variation 
of  the  nominal  axial  velocity  V^^  for  r/Rp  =  .7.  First,  we  consider  the 
results  without  V^.  expected,  based  on  r(9p)  and  nearly 

symmetric  with  respect  to  its  maximum  value  at  9p  =  -  24®.  K^^(9p),  which 

by  definition  is  almost  proportional  to  dP/dt,  is  nearly  antisymmetric  with 
respect  to  9p  =  0.  K.j.^  and  combine  such  that  the  total  thrust  K-p  is 

also  nearly  symmetric  with  respect  to  its  peak  value  at  9p  =  8®,  and  fur¬ 

thermore,  can  be  directly  correlated  as  a  quasi-steady  response  to 
Although  the  results  including  V  are  similar,  there  is  a  remarkable  dif- 
ference  such  that  K,p  and  are  no  longer  symmetric,  but  indicate  signifi¬ 

cantly  larger  values  on  the  starboard  than  on  the  port  side  of  the  hull. 
This  important  influence  of  V  on  the  propeller  loading  will  be  discussed 
further  below. 

The  experiments  described  above  Involved  steady  mean-flow  measure¬ 
ments,  i.e.,  averages  were  taken  over  a  sufficient  length  of  time  to  remove 
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the  unsteadiness  of  both  the  turbulence  and  the  propeller  to  define  tine- 
mean  values.  Therefore,  to  relate  the  propeller  performance  to  the  mean- 
flow  measurements,  of  Importance  is  not  the  unsteady  blade  thrust  and  tor¬ 
que,  but  rather  the  time-averaged  angular  variation  of  thrust  and  torque. 
As  mentioned  above,  the  propeller-performance  program  was  modified  to  pro¬ 
vide  the  latter.  The  usual  definitions  for  unsteady  blade  thrust  and  tor¬ 
que  are 
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and  5  “  9  ~  9p  and  9p  =  2irnt  (see  List  of  Symbols  for  definitions).  On  the 
other  hand,  the  definitions  for  the  time-averaged  angular  variation  of 
thrust  and  torque  are 

lLp(9)  =  f  f^(r,e)rdc  (IV-6) 

K^(9)  “  f  (r,0)r^dr  (lV-7) 

where 
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The  steady  thrust  and  torque  are  obtained  equivalently  by  either  time- 
averaging  (i.e.,  averaging  over  9p)  (IV-3)  and  (IV-4),  or  circumferenti¬ 
ally-averaging  (i.e.  ,  averaging  over  9)  (IV-6)  and  (lV-7).  The  calculated 

values  based  on  the  nominal  inflow  are  compared  with  the  present  measure¬ 
ments  in  Table  5.  It  is  seen  that  the  calculated  values  overpredict  the 
thrust  and  torque  by  about  13%  and  8%,  respectively.  Also  shown  in  Table  5 
are  the  calculated  steady  side  forces  and  moments  as  well  as  the  three 
components  of  the  unsteady  first  blade  harmonic  forces  and  moments, 
although  none  of  these  quantities  were  measured. 

Figure  26  shows  a  comparison  of  the  unsteady  and  time-averaged  angular 
variation  of  blade  thrust  and  torque.  For  comparison  purposes,  K^(9)  and 
Kq(9)i  which  represent  total  thrust  and  torque,  are  converted  to  per  blade 
values  by  the  multiplicative  factor  2tt/N,  where  N  is  the  number  of  blades. 
It  is  seen  that  both  the  unsteady  and  time-averaged  angular  varying  thrust 
and  torque  are  quite  similar.  The  main  difference  is  that  the  time-aver¬ 
aged  angular  varying  values  are  shifted  towards  the  1/4  chord  (forward  of 
the  blade  generator  line,  see  figure  4b)  relative  to  the  unsteady  values. 
We  also  see  significantly  larger  loading  on  the  starboard  than  the  port 
side  of  the  hull  for  reasons  mentioned  earlier.  To  display  further  the 
characteristics  of  the  propeller  loading,  figure  27a, b  shows  the  time- 
averaged  angular  and  radially  varying  thrust  and  torque  contours,  i.e., 
f  (r,9)  and  f  (r,9),  defined  in  (IV-8)  and  (IV-9),  respectively.  Also 

X  9 

shown  in  parentheses  on  figure  27a  are  the  corresponding  thrust-loading 
2  2  —  2 

values  (i.e.,  2n  Dp  ).  Clearly  evident  are  the  large  values  of  thrust 

for  0  <  9  <  90°  and  near  9  =  -150°  (figure  27a).  Also,  note  that,  in  the 
lower  part  of  the  propeller  disk  near  the  hub,  the  thrust  is  generally 
larger  on  the  starboard  than  the  port  side  of  the  hull,  in  this  case,  due 
to  the  influence  of  the  weak  longitudinal  vortex.  The  torque  contours  are 
similar  to  those  for  thrust,  but  with  greater  symmetry  and  more  gradual 
variations  (figure  27b). 

Lastly,  regarding  the  propeller  performance,  consideration  is  given  to 
the  results  for  the  effective  inflow.  As  mentioned  above,  the  effective 
inflow  was  calculated  using  the  method  of  Toda  et  al.  (1987)  (see  Appendix 
I).  In  this  approach,  the  method  of  Huang  and  associates  (1976,  1980)  is 
applied  two-dimensionally  in  x-y  planes  (i.e.  ,  for  z  =  constant)  to  obtain 
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the  lateral  propeller-induced  flow  contraction.  The  vertical  contraction 
is  determined  based  on  that  for  the  propeller  operating  in  open  water  with 
the  same  advance  coefficient  as  the  model  (or  ship).  For  the  present 
application,  only  the  axial  component  of  the  effective  inflow  was  calcu¬ 
lated,  since  for  fine  hull  forms,  the  effects  of  the  propeller  on  the 
crossplane  components  for  the  region  upstream  of  the  propeller  are  quite 
small  (see  figures  21f-h,  also  Kasahara  (1985)).  As  will  be  shown  below, 
the  primary  differences  between  the  nominal  and  effective  inflows  are  that 
the  latter  has  increased  velocity  and  uniformity. 

Included  in  figure  2^  are  the  circumferential-average  axial  velocity 

(figure  24a)  and  the  spanwii,  circulation  distribution  r  (figure  24b)  for 

the  effective  inflow.  Although  the  shape  of  the  distribution  for  the 

effective  Inflow  is  similar  to  that  for  the  nominal  inflow,  there  is  a 

slight  shift  towards  the  uniform-inflow  distribution.  With  the  effective 

inflow,  r  is  smaller  than  that  for  either  the  nominal  or  uniform  inflow 
max 

and  occurs  at  r/Rp  -  .64.  These  differences  between  the  results  for  the 
nominal  and  effective  Inflows  are  due  to  the  increased  uniformity  of  the 
latter  resulting  in  smaller  section  angle  of  attack,  variations.  This 
effect  is  clearly  demonstrated  by  the  circumferential-average  axial  veloc¬ 
ity  (figure  24a),  which  also  shows  that  the  largest  increases  in  velocity 
for  the  effective  inflow  are  for  the  inner  radii. 

Figure  28  shows  a  comparison  of  the  unsteady  and  time-averaged  angular 
variation  (converted  to  per-blade  values)  of  blade  thrust  and  torque  for 
the  effective  inflow.  Also  shown  is  the  variation  of  the  effective  axial 
velocity  at  r/Rp  *  .7  as  well  as  that  for  the  nominal  inflow,  reproduced 
here  from  figure  25  for  the  purpose  of  direct  comparison.  Quite  apparent 
is  the  already  pointed  out  increased  velocity  and  uniformity  of  the  effec¬ 
tive  as  compared  to  the  nominal  inflow.  The  velocity  increase  is  nearly 
the  same  for  all  angular  positions  (0  or  9p).  A  comparison  of  figures  26 
and  28  indicates  that  the  results  for  the  effective  inflow  are  similar  and 
consistent  with  those  described  earlier  for  the  nominal  inflow,  but  with 
decreased  values  and  differences  between  the  starboard  and  port  sides  of 
the  hull  for  reasons  mentioned  earlier.  Also  included  in  Table  5  are  the 
calculated  values  of  steady  and  first  blade  harmonic  forces  and  moments 
based  on  the  effective  inflow.  It  is  seen  that,  in  this  case,  the  calcula- 
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tions  for  the  steady  thrust  and  torque  are  in  excellent  agreement  with  the 
measurements. 

Finally,  for  the  effective  Inflow,  figure  29  shows  the  time-averaged 
angular  and  radially  varying  thrust  f^(r,9)  (figure  29a)  and  torque 
f.(r,9)  (figure  29b)  contours  using  the  same  format  as  figure  27.  A 
comparison  of  figures  27  and  29  indicates  that,  here  again,  the  results  for 
the  effective  inflow  are  similar  and  consistent  with  those  described 
earlier  for  the  nominal  inflow,  but  with  decreased  values,  especially  near 
the  hull  centerplane,  for  reasons  mentioned  earlier.  For  the  present 
application,  it  appears  that  the  effective  and  nominal  inflows  are  qualita¬ 
tively,  but  not  quantitatively,  similar. 

Keeping  in  mind  the  above  discussion  of  propeller  performance,  con¬ 
sideration  is  now  given  to  propeller-hull  interaction.  First,  we  consider 
the  interaction  for  the  region  upstream  of  the  propeller  plane.  In  order 
to  aid  in  explaining  the  propeller  effects  pointed  out  earlier,  figure  30 
shows  contours  of  the  difference  in  axial  velocity  Au  between  the  with- 
and  without-prope^ler  conditions  for  x  =  .975  and  .98125.  Also  shown  is 
Au  vs.  y  at  various  horizontal  planes.  The  circumf erential-averaged 
propeller-induced  axial  velocity  for  both  the  nominal  and  effective 

inflows  is  included  for  comparison.  Unfortunately,  our  version  of  the 
propeller-performance  program  only  allows  for  the  evaluation  of  steady 
(i.e.,  circumferentially  averaged)  field-point  velocities.  However,  for 
the  region  upstream  of  the  propeller  plane,  unsteady  effects  on  the  field- 
point  velocities  are  negligible.  It  is  apparent  that  the  low-momentum 
fluid  near  the  hull  and  in  the  longitudinal-vortex  core  undergo  the  largest 
axial  acceleration.  Within  the  propeller  disk,  Au  is  considerably  larger 
than  Vp^  for  both  inflows,  which  indicates  that  the  increase  in  velocity  is 
not  simply  due  to  the  superposition  of  the  propeller-induced  flow  and  the 
hull  boundary  layer,  but  the  result  of  the  contraction  and  axial  acceler¬ 
ation  of  the  hull  boundary  layer  due  to  the  action  of  the  propeller.  The 
magnitude  of  the  axial  acceleration  depends  both  on  hull  proximity  and  on 
the  distribution  of  propeller  loading,  i.e.,  radial  location.  The  flow  is 
very  nearly  symmetric  with  respect  to  the  hull  centerplane.  Consistent 
with  the  previous  discussion,  the  Vpx  values  for  both  the  nominal  and 
effective  inflow  are  quite  similar,  but  with  the  latter  slightly  smaller 
due  to  its  reduced  loading. 
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Next,  we  consider  the  interaction  in  the  propeller  slipstream.  As 
will  now  be  discussed,  most  of  the  flow-field  features  pointed  out  earlier 
can  be  explained  as  a  direct  consequence  of  the  time-averaged  angular  and 
radially  varying  propeller  loading.  A  comparison  of  the  loading  contours 
(figures  27  and  29)  with  both  the  axial  velocity  (figure  21j,k)  and  pres¬ 
sure  (figure  23j )  contours  in  the  very  near  wake  indicates  that  all  three 
have  a  very  similar  pattern,  i.e. ,  increased  values  on  the  starboard  side 
of  the  hull  for  0°  <  9p  <  90°  and  on  the  port  side  for  ©p  =  -150°  in  com¬ 
parison  to  the  other  regions  within  the  propeller  disk.  Note  that  the 
increased  values  of  pressure  on  the  starboard  as  compared  to  the  port  side 
of  the  hull  are,  no  doubt,  the  cause  of  the  drift  of  the  wake  centerplane 
off  the  hull  centerplane  towards  the  port  side.  Subsequently,  in  the  near 
wake,  the  axial  velocity  (figures  2l£,m)  and  pressure  (figures  23k-o) 
contours  become  more  uniform  due  to  diffusion.  The  rate  of  recovery  of  the 
pressure  is  considerably  larger  than  that  of  the  axial  velocity.  Thus,  it 
appears  that  the  propeller  loading  resulting  from  its  operation  in  nonuni¬ 
form  Inflow  is  responsible  for  the  nature  of  the  complicated  flow  within 
the  propeller  slipstream.  This  effect  is  further  displayed  in  figure  31 
which  shows  the  difference  in  axial  velocity  contours  between  the  with-  and 
without-propeller  conditions  Au  for  x  =  1.  Also  shown  is  Au  vs.  y  for 
various  horizontal  planes,  including  the  time-averaged  angular  and  radially 
varying  thrust  f^  and  propeller-induced  axial  velocity  for  both  the 
nominal  and  effective  inflows  for  comparison.  Note  that  f^  and  Vp^  are 
evaluated  at  the  propeller  plane.  Two  methods  were  used  to  calculate 
Vp^.  The  first  method  is  the  propeller-performance  program  which,  as 
already  mentioned,  provides  circumferentially-averaged  values.  The  second 
method  is  based  on  the  momentum  theorem,  i.e. 

V  ,  ■  - - - 

'^Px  =  (-1  +  /I  +  (IV-10) 

The  latter  method,  although  quite  approximate,  includes  unsteady  effects, 
which  are  known  to  be  important  in  the  propeller  slipstream.  A  comparison 
of  figures  31  and  27  or  29  indicates  that  the  Au  contours  show  even  greater 
similarity  with  the  propeller  loading  than  that  noted  above  for  the  axial 
velocity  and  pressure  contours.  For  z  “  .02  and  .045  the  shape  of  the  Au 
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For  z 


vs.  y  curves  are  nearly  identical  to  those  of 


.033,  the  dif¬ 


ferences  between  Au  and  f^  are  due  to  the  influence  of  the  hub  and  its 
vortex.  Just  as  was  the  case  for  the  region  upstream  of  the  propeller 
plane,  within  the  propeller  disk,  Au  is  considerably  larger  than  Vp^,  again 
indicating  a  significant  acceleration  of  the  hull  boundary  and  wake  due  to 
the  action  of  the  propeller.  However,  in  this  case,  only  a  qualitative 
assessment  is  possible  due  to  the  approximations  involved  in  calculating 


'Px 


and  due  to  the  fact  that  the  values  shown  are  for  the  propeller  plane 


and  not  x  =  1.  Note  that  the  interaction  is  largest  near  the  hull  center- 
plane.  Also,  for  z  =  .02  the  boundary  layer  is  thick  such  that  the  f^  and 
Au  curves  are  fairly  broad,  whereas  for  z  =  .045  the  boundary  layer  is 
thin  and  the  f^  and  Au  curves  are  quite  sharply  peaked.  Here  again,  the 
Vpx  values  for  both  the  nominal  and  effective  inflow  are  quite  similar. 
The  differences  are  consistent  with  the  previous  discussions. 


V.  CONCLUDING  REMARKS 

Detailed  experimental  information  has  been  presented  which  documents 
the  interaction  between  a  propeller  and  the  hull  boundary  layer  and  wake 
for  a  Series  60  Cg  =  .6  model  ship.  The  experimental  equipment  and  proce¬ 
dures  have  been  described,  and  the  results  discussed  to  point  out  the  es¬ 
sential  differences  between  the  flows  with  and  without  propeller.  The 
results  have  been  analyzed  to  assess  the  nature  of  the  interaction  between 
the  propeller  and  the  hull  boundary  layer  and  wake.  To  this  end,  use  is 
made  of  a  propeller-performance  program  with  both  nominal  and  effective 
inflows.  The  interaction  for  the  region  upstream  of  the  propeller  is 
mainly  the  result  of  the  contraction  and  axial  acceleration  of  the  hull 
boundary  layer  due  to  the  action  of  the  propeller.  Thus,  for  the  present 
application,  a  simple  prediction  method  for  the  effective  Inflow  was  shown 
to  be  adequate.  However,  for  more  complex  geometries  (e.g.,  full  hull 
forms,  appended  bodies,  etc.)  it  is  expected  that  more  comprehensive 
methods  will  be  required.  The  interaction  in  the  propeller  slipstream  is 
quite  complicated,  i.e.,  the  flow  is  completely  altered  due  to  the  action 
of  the  propeller,  indicating  characteristics  which  are  similar  to  a  swirl¬ 
ing  jet.  However,  most  features  can  be  explained  as  a  direct  consequence 
of  the  propeller  loading  resulting  from  its  operation  with  a  three- 
dimensional  nonuniform  Inflow. 
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Although  we  have  been  able  to  qualitatively  explain  most  of  the  exper¬ 
imental  results  solely  with  reference  to  the  propeller  loading  predicted 
from  a  propeller-performance  program  with  the  nominal  and  effective 
inflows,  it  should  be  recognized  that  quantitative  prediction  of  the  com¬ 
plete  flow  field  requires  the  use  of  advanced  viscous-flow  methods,  includ¬ 
ing  the  effects  of  the  propeller  (Stern  et  al.,  1988a).  In  fact,  one  of 
the  motivations  of  the  present  study  was  to  provide  documentation  of  both 
the  test  conditions  and  results  in  sufficient  detail  to  be  useful  as  a  test 
case  for  validating  computational  methods.  The  data  discussed  here  are 
available  on  magnetic  tape  upon  request  from  IIHR.  A  comparison  of  the 
present  experimental  results  with  the  predictions  of  the  method  of  Stern  et 
al.  (1988a)  is  in  progress  and  will  be  reported  in  the  near  future. 

Finally,  with  regard  to  the  direction  of  future  work,  additional 
experiments  in  the  intermediate-  and  far-wake  regions  are  desirable  to 
determine  the  evolution  and  recovery  of  the  wake  as  well  as  experiments  at 
higher  Fr  to  determine  free-surface  effects.  Also,  of  Interest  are  exper¬ 
iments  for  other  geometries,  e.g.  ,  full  hull  forms,  high-speed  ships,  and 
appended  bodies. 
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Table  1 

Principal  dimensions  of  ship  model  and  propeller 


model  propeller 


Length 

4.0m 

Diameter  (=  2Rp) 

145.64mm 

Breadth 

0.5333m 

Pitch  ratio  (constant  pitch) 

1.0310 

Draft 

0.2133m 

Boss  ratio 

0.2 

Cb 

0.6 

Number  of  blades 

5 

Cm 

.977 

Direction  of  turning 

right 

m  *  (A^/Aj) 

.0033 

Expanded  area  ratio 

0.7451 

Zq  (shaft  center) 

129.05mm 

Projected  area  ratio 

0.6249 

Wetted  surface  area 

2.7189m2 

Thickness  ratio  at  0.7Rp 

3.992% 

Thickness  distribution 

n  =  25 

Blade  section 

MAU 

Rake 

6° 

Table  2 

Measurement  locations  for  the  with-  and  without-propeller 
conditions  (port  and  starboard) 


X  x/Rp  Surface-  Mean-  WT  Re 

Pressure  Velocity  ®C  x 

Distribution  and  10“^ 

Pressure 
Field 


.5 

-26.79 

w/ 0 

19.5 

3.94 

.6 

-21.29 

w/  o 

19.5 

3.94 

.7 

-15.8 

w/  0 

19.5 

3.94 

.8 

-10.3 

w,w/o 

w/ o 

10.6 

3.11 

.85 

-7.55 

w,w/o 

.9 

-4.81 

w,w/  o 

w,w/o 

10.2 

3.08 

.925 

-3.43 

w,w/o 

.95 

-2.06 

w,w/  o 

w,w/o 

8.0 

2.89 

.9625 

-1.37 

w,w/o 

.975 

-.69 

w,w/  o 

w,  w/o 

11.6 

3.20 

.981 

-.36 

w,w/ o 

.98125 

-.34 

w,w/o 

8.0 

2.89 

Propeller- 

.9875 

.0 

w,w/o 

w/o 

3.0 

2.89 

Plane 

1. 

.69 

w,w/o 

16.0 

3.60 

1.01875 

1.72 

w,  w/o 

10.7 

3.12 

1.05 

3.43 

w,w/o 

15.4 

3.55 

1.1 

6.18 

w,  w/o 

14.9 

3.50 

Resistance  test  (IKb^C) 

Self-propulsion  test  (ll.7‘’C) 

Propeller-open  water  test  (11.5'’C) 

Surface-Pressure  Distribution  (10.2‘’C,  Re  »  3.08  x  10  ) 
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Table  3 
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i) 
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JTJ 


OJ 
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♦- 

oo 
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(N 

O  w-N 
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3  4) 
II 

^  C 

(X)  4) 

OLU 


/XJ_Opfl3 

w-  r3 
0(Diy) 


ese 

^  sO  00 
<N 

ON  CM 


ooc  ^ 

c  fo 

1)  4;  4.^ 
^cOO 


(U 

c 


rtj 


•Offsets  are  in  meters  for  full  scale 


Table  <( 

Propeller  ollsets 
Propeller  geometry 


r/R 

0.2 

0.3 

0.4 

0.5 

0.6 

0.66 

0.7 

0.8 

0.9 

1.0 

Chord  Lenj-th  ^ 

^0.66 

From  Generator 

27.95 

33.45 

38.76 

43.54 

47.96 

49.74 

51.33 

52.39 

48.49 

17.29 

Line  to  Trailing 
Edge 

(Chord  length  at  .66R 

Co.66  = 

From  Generator 
Line  to  Leading 
Edge 

38.58 

44.25 

48.32 

50.80 

51.15 

50.26 

48.31 

40.53 

25.13 

Chord  Length 

66.54 

77.70 

87.08 

94.34 

99.11 

100.00 

99.64 

92.92 

73.62 

- 

Max.  Thickness 

Di aneter 

4.06 

3.59 

3.12 

2.65 

2.18 

1.90 

1.71 

\.2U 

0.77 

0.30  B.T.F.  =  0.050 

Maximum  thickness 

Position  from  L.E. 

32.0 

32.0 

32.0 

32.5 

34.9 

37.9 

40.2 

uyi* 

48.9 

- 

Wing  section  oHsets 

<:  distance  (rom  L.E.  (%  of  chord) 

I ;  %  ot  maximum  thickness;  =  upper  surface;  yj  =  lower  surface 


r/R  X 
0.20 

X 

0 

35.00 

2.00 

51.85 

24.25 

4.00 

59.75 

19.05 

6.00 

66.15 

15.00 

10.00 

76.05 

10.00 

15.00 

85.25 

5.40 

20.00 

92.20 

2.35 

30.00 

99.80 

32.00 

100.00 

40.00 

97.75 

50.00 

89.95 

60.00 

78.15 

70.00 

63.15 

80.00 

45.25 

90.00 

25.30 

95.00 

15.00 

100.00 

4.50 

X 

0 

2.00 

4.00 

6.00 

10.00 

15.00 

20.00 

30.00 

32.00 

40.00 

50.00 

60.00 

70.00 

80.00 

90.00 

95.00 

100.00 

0.  30 

?? 

35.00 

51.85 

24.25 

59.75 

19.05 

66.15 

15.00 

76.05 

10.00 

85.25 

5.40 

92.20 

2.35 

99.80 

100.00 

97.75 

89.95 

78.15 

63.15 

45.25 

25.30 

15.00 

4.50 

0,40 

X 

0 

2.00 

4.00 

6.00 

10.00 

15.00 

20.00 

30.00 

32.00 

40.00 

50.00 

60.00 

70.00 

80.00 

90.00 

95.00 

100.00 

35.00 

51.85 

24.25 

59.75 

19.05 

66.15 

15.00 

76.05 

10.00 

85.25 

5.40 

92.20 

2.35 

99.80 

100.00 

97.75 

89.95 

78.15 

63.15 

45.25 

25.30 

15.00 

0.50 

0.30 

X 

0 

2.03 

4.06 

6.09 

10.16 

15.23 

20.31 

30.47 

32.50 

40.44 

50.37 

60.29 

70.22 

80.15 

90.07 

95.04 

100.00 

Jr 

35.00 

24.25 

51.85 

19.05 

59.75 

15.00 

66.15 

10.00 

76.05 

5.40 

85.23 

2.35 

92.20 

99.80 

100.00 

97.75 

89.95 

78.15 

63.15 

45.25 

25.30 

15.00 

4.50 

0,60 

X 

0 

2.18 

4.36 

6.54 

10.91 

16.36 

21.81 

32.72 

34.90 

42.56 

52.13 

61.70 

71.28 

80.85 

90.43 

95.21 

100.00 

Jr 

34.00 

49.60 

23.60 

58.00 

18.10 

64.75 

14.23 

75.20 

9.45 

84.80 

5.00 

91.80 

2.23 

99.80 

100.00 

97.75 

89.95 

78.15 

63.15 

45.25 

25.30 

15.00 

a. 50 

0.70 

X 

0 

2.51 

5.03 

7.54 

12.56 

18.84 

25.12 

37.69 

34.20 

47.23 

56.03 

64.82 

73.62 

82.41 

91.21 

95.60 

100.00 

Jr 

30.00 

42.90 

20.50 

32.20 

15.45 

59.90 

11.95 

71.65 

7.70 

82.35 

4.10 

90.60 

1.73 

99.80 

100.00 

97.75 

89.95 

78.15 

63.15 

45.25 

25.30 

15.00 

0.50 

0.80 

X 

0 

2.51 

5.68 

8.51 

14.19 

21.28 

28.38 

42.56 

45.40 

51.82 

59.85 

67.89 

75.91 

83.94 

91.97 

95.99 

100.00 

Jr 

21.00 

32.45 

14.00 

41.70 

10.45 

50.10 

8.05 

64.60 

5.05 

78.45 

2.70 

88.90 

1.15 

99.80 

100.00 

97.75 

89.95 

78.15 

63.15 

45.23 

25.30 

15.00 

0.90 

X 

0 

3.06 

6.11 

9.17 

15.28 

22.92 

30.56 

45.85 

48.90 

54.91 

62.42 

69.94 

77.46 

84.97 

92.49 

96.24 

100.00 

Jr 

8.30 

21.10 

4.00 

31.50 

2.70 

40.90 

2.05 

57.45 

1.20 

74.70 

0.70 

87.45 

0.30 

99.70 

100.00 

98.65 

92.73 

83.00 

69.35 

51.85 

30.80 

I9.«i0 

6.S4 

0.95 

X 

0 

3.13 

6.25 

9.38 

15.63 

23.44 

31.25 

46.87 

50.00 

55.88 

63.23 

70.59 

77.94 

85.30 

92.65 

96.32 

100.00 

Jr 

6.00 

19.65 

30.00 

39.60 

56.75 

74.30 

87.30 

99.65 

100.00 

99.00 

93.85 

84.65 

71.65 

34.30 

33.50 

21.50 

s.oo 
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Table  5 

Steady  and  unsteady  first  blade  harmonic  propeller 
forces  and  moments 

steady 


nominal  effective 

inflow  inflow  experiment 


Forces 

Thrust  (K.j>) 

Horizontal  (+  port) 
Vertical  (+  upward) 

.265 

.00244 

.00982 

.233 

.00403 

.00962 

.234 

Torque  (Kg) 

.0444 

.0403 

.0411 

Moments 

Horizontal  (+  port) 

.00793 

.00886 

Vertical  (+  upward) 

.00673 

.00657 

unsteady 

first  blade  harmonic 

nominal 

effective 

inflow 

inflow 

Thrust 

.00665 

.00670 

Forces 

Horizontal  (+  port) 

.00858 

.00837 

Vertical  (+  upward) 

.00440 

.00463 

Torque 

.00080 

.00090 

Moments 

Horizontal  (+  port) 

.00771 

.00693 

Vertical  (+  upward) 

.00205 

.00214 

Forces  are  nondlmensionalized  using  pn  Dp 

2  5 

Moments  are  nondlmensionalized  using  pn  Dp 
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es 


Figure  2.  Osaka  University,  Department  of  Naval  Architecture,  towing  tank. 


0  8  x/L«0.9 
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(b)  profile  plan  (longitudinal  section  through  centerplane) 
Figure  3.  Series  60  Cg  =  .6  lines  drawing. 
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(b)  coordinate  system 

Propeller  drawing  and  coordinate  system. 
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Figure  8.  Block  diagrams  for  mean-velocity  and 
pressure  field  measurements. 
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m/5 


Figure  16.  Self-propulsion  factors:  (a)  wake  fraction  (1  w) ;  (b)  thrust 
deduction  (1-t);  and  (c)  relative  rotative  efficiency 
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FiRure  17.  Wave  profiles. 
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velocity  crossplane  vectors  and  axial  velocity  contours  (continued) 
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Figure  21.  (continued) 
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Figure  22.  (continued) 
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Figure  30.  Difference  in  axial  velocity  contours  between  the  with 
and  without-propeller  conditions  au  for  the  region  up¬ 
stream  of  the  propeller. 


Appendix  I.  English  Translation  of  Toda  et  al.  (1987) 

"A  Simple  Prediction  Method  for  Effective  Wake  Distribution” 

by 


Y.  Toda^,  Y.  Kasahara^,  and  I.  Tanaka^ 


Abstract 

The  interaction  between  a  propeller  and  a  flow  field  around  a  full  form 
ship  is  investigated.  A  method  is  presented  for  the  prediction  of  effective 
velocity  field  from  measured  nominal  velocity  field.  The  prediction  is  made 
two-dimensionally ,  with  a  three-dimensional  correction,  by  use  of  a  modified 
form  of  Huang  et  al.'s  method  under  the  assumption  that  the  effects  of  the 
vorticity,  whose  axis  is  vertical,  is  dominantly  responsible  for  the  distor¬ 
tion  of  velocity  distribution  due  to  propeller  suction.  The  agreement  between 
measured  and  predicted  total  velocity  distribution  with  propeller  is  compar¬ 
atively  good.  Characteristics  of  the  effective  velocity  field  for  a  full  form 
ship  are  discussed. 

1.  Introduction 

In  the  flow  field  around  a  ship  stern  where  a  propeller  is  operating, 
there  exists  an  interaction  between  the  viscous  flow  field  and  the  propel¬ 
ler.  The  flow  field  where  the  propeller  is  operating  cannot  be  understood  as 
a  superposition  of  propeller  induced  velocity  field  and  nominal  velocity 
field,  which  is  the  inflow  velocity  distribution  at  the  propeller  plane  in  the 
absence  of  the  propeller.  This  is  why  the  strength  and  the  location  of  the 
vortices  in  boundary  layer  are  changed  due  to  propeller  suction.  At  present, 
the  calculation  methods  for  propeller  performance,  propeller  bearing  forces, 
surface  forces,  and  so  on  are  considered  to  be  mostly  those  taking  the  veloc¬ 
ity  distribution  at  the  propeller  plane  as  the  input  and  satisfying  the  bound¬ 
ary  condition  by  adding  induced  velocity.  Therefore,  it  is  considered  that  as 
the  input  for  propeller  calculation,  effective  velocity  distribution  should  be 
used  instead  of  nominal  velocity  distribution.  Effective  velocity  field  is 
considered  to  be  the  flow  field  induced  by  the  distorted  vorticity  distribu¬ 
tion  in  nominal  flow  field  due  to  propeller  suction. 

Recently  various  works  on  effective  velocity  have  been  carried  out.  In 
particular,  as  for  axisymmetric  flow,  a  lot  of  works  have  been  done  (for 
example,  [  1  ]  ,  [ 2 ] , [ 3 ] , [ 4 ] )  and  investigations  utilizing  the  Reynolds  equations 
have  also  been  performed  (for  example,  [5], [6]).  As  for  a  simple  three  dimen¬ 
sional  hull  iorm,  the  authors  [7],  Sato  et  al.  [8],  and  others  calculated  the 
effective  velocity  distribution  by  using  an  integral  method  for  calculating 
the  boundary  layer.  As  for  practical  hull  forms,  Breslin  et  al.  [9],  Chen  et 


1  Osaka  University 

2  Tsu  Research  Laboratories,  Nippon  Kokan,  K.K. 
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al.  [10],  Dyne  [11],  van  Gent  [12],  and  others  have  shown  methods  for  predict¬ 
ing  the  effective  wake  distribution  from  the  nominal  velocity  distribution  by 
methods  extended  from  the  axisymmetric  case. 

This  paper  describes  a  simple  method  for  calculating  the  effective  veloc¬ 
ity  field  from  the  nominal  velocity  field  by  use  of  a  modified  form  of  Huang 
et  al.'s  method  [1].  The  prediction  is  made  two-dinensionally ,  with  a  three- 
dimensional  correction  under  the  assumption  that  the  effects  of  the  change  of 
vorticity,  whose  axis  is  vertical,  is  dominantly  responsible  for  the  distor¬ 
tion  of  velocity  distribution  due  to  propeller  suction.  The  calculation  is 
made  for  a  full  form  ship.  The  calculated  velocity  distribution  with  propel¬ 
ler  in  the  region  ahead  of  the  propeller  is  compared  with  measured  results. 
The  effective  velocity  distribution  at  the  propeller  plane  is  calculated  and 
the  difference  between  the  nominal  and  effective  velocity  distributions  are 
investigated.  The  propeller  thrust  is  also  calculated  using  the  effective 
velocity,  as  well  as  using  the  nominal  velocity  and  they  are  compared  with  the 
measured  thrust. 


2.  A  prediction  method  for  effective  velocity  distribution 

As  shown  in  Fig.  1,  we  adopt  a  right-handed  Cartesian  coordinate  sys¬ 
tem.  The  x-axis  coincides  with  the  propeller  axis  of  revolution  and  the 
direction  of  uniform  flow,  the  y-axis  is  in  the  starboard  direction,  and  the 
z-axis  is  in  the  upward  direction.  Let  u,  v  and  w  be  the  nondimensional 
velocity  components  ( nondiraensionalized  by  the  uniform  flow)  in  the  direction 
of  the  X,  y  and  z  axes,  respectively. 

In  this  paper,  it  is  assumed  that  the  effect  of  the  change  of  vorticity 
component  ,  whose  axis  is  vertical,  is  dominantly  responsible  for  the  dis¬ 
tortion  of  the  velocity  distribution  due  to  propeller  suction  (from  the  mea¬ 
sured  flow  fields  around  a  full  form  ship  both  with  and  without  propeller 
[13]).  Therefore,  the  change  of  oj^  is  mainly  considerad  two-dimensionally  for 
the  velocity  distribution  on  a  line  where  x  and  z  are  constant  and,  as  for  the 
three-dimensional  effects,  the  flow  contraction  in  the  z  direction  in  irrota- 
tional  flow  is  taken  into  account  in  a  similar  manner  as  in  a  previous  paper 
[7]. 


First,  it  is  simply  shown  how  Huang  et  al.'s  method  [1]  for  an  axisym¬ 
metric  case  transforms  for  a  two-dimensional  case.  It  is  assumed  that  fluid 

is  inviscid  but  rotational  similarly  to  Huang  et  al.  In  the  case  of  axisym¬ 
metric  flow,  Wg/r  is  constant  along  streamlines,  while  in  th^  case  of  two- 
dimensional  flov/,  0)  is  constant  along  streamlines.  The  nominal  velocity 
without  propeller  and  the  total  velocity  with  propeller  are  denoted  by  u  , 
and  u  ,  respectively.  The  total  velocity  is  divided  into  the  effective  veloc¬ 
ity  u^  and  the  induced  velocity  u^^Cu  =u  +u.;  suffixes  are  used  similarly  for  v 
and  w^.  Also,  as  for  the  coordinate^y,®sulf ixes  p  and  n  denote  the  conditions 
with  and  without  propeller,  respectively.  Since  m  is  constant  on  a  stream¬ 
line  (see  Fig.  2),  ^ 
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is  obcained.  The  induced  velocity  field  is  assumed  to  be  irrotational.  The 
usual  boundary-layer  approximation  9v/3x  «  9u/ 3y  is  assumed  to  be  valid  for 
both  the  nominal  velocity  and  the  effective  velocity.  Then  Eq.  (1)  reduces  to 


The  massflow  condition  is  written  as 


u  dy 
n  n 


u  dy 
P  F 


=  (u  +  u  )dy 
e  i  t 


From  Eq.  (1)  and  Eq.  (3) 


u  du  =  (u  +  u . )dy  (4) 

n  n  e  i  P 

is  obtained.  Eqs.  (3)  and  (4)  are  the  governing  equations  for  the  interaction 
between  the  propeller  and  the  stern  flow  field.  If  the  distributions  u^  and 
u-  are  given,  the  distributions  of  u  and  u  can  be  determined.  (The  norainal 

J.  ®  ’  P  0 

velocity  at  the  surface  is  taken  to  be  the  extrapolated  value  of  the  nominal 
velocities  from  near  the  surface.) 

Next,  the  finite  difference  form  of  Eqs.  (3)  and  (4)  are  discussed.  A.s 
shown  in  Fig.  3,  if  the  distribution  of  u  on  a  cross  section  of  x-constant  is 
given  by  the  measured  one  or  others,  tKe  values  at  the  n  discrete  points 

^nl’^n2’  ^n3’  ’  *  *  ’  ^no  denoted  by  »'^n2 ''^n3  ’  *  *  ’  »^nn’  >^espectively . 

Also  for  otner  quantities,  the  values  at  discrete  points  are  denoted  similarly 

by  second  suffixes.  Then  Eqs.  (3)  and  (4)  can  be  written  as 


( u  , +u  .,,)(u  u  .)  =  (u  .+  U..+  u  u...,)(u  u  .) 

'  ni  ni+1  ni+1  ni  '  ei  ii  ei+1  ii+l^^  ei+1  ei' 

(u  .+  u  ._i.,)(y  V  .)  =  (u  .+  u.  .+  u  u..,,)(y  y  .) 

ni  ni+l  ni+1  ni  ei  ii  ei+1  ii+l  "^pi+l  pi 


>(^ei+r  '^ei^ 

(5) 

.^^^pi+r  ^pi^ 

(6) 

Since  the  boundary-layer  approximations  are  assumed  to  be  valid  for  the  nom¬ 
inal  and  the  effective  velocity,  the  vorticity  within  the  boundary  layer  does 
not  affect  the  velocity  field  outside  the  boundary  layer  and  it  is  considered 
that  Ug  =  Uj^  outside  the  boundary  layer.  Therefore,  if  y^^  is  taken  to  be 
outside  the  boundary  layer,  the  outside  edge  condition  is  written  as 


u  =  u 
en  nn 
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Or,  if  the  effect  of  the  irrotational  flow  is  taken  into  account,  the  outside 
edge  condition  is  written  as 


u  (y  )=u  (v  )-(u  (v  )-u  (y  )) 
en^-'pn'^  nn^-’^nn'^  pot  nn  pot'-'pn^^ 


(8) 


where,  u  is  the  value  obtained  by  potential  calculation.  The  body  surface 
is  a  streamline.  Therefore  the  inside  edge  condition  is  written  as 


^nl  =  ypl  ° 


(9) 


where,  y^  is  the  y-coordinate  of  the  body  surface  or  the  center  plane.  By 
simply  changing,  Eq.  (5)  reduces  to 


2  2  ,  11  11+ K2  , 

u  u  .=  (u  +  - = - )  -  (u  + 

ni+1  ni  ei+1  2  ei 


u  .  .+  u .  .  ,  ,  « 
11  ii+l,^2 

'  n  f 


(10) 


u  . 
ei 


u .  .  +  u. .  .  ,  „  « 

//  ,  11-  11+1^2  .  t  2 

+ - - - )  +  (u  . 

ei+1  2  ni 


2  ,  ^ii+1 

-  - 


+  u. 


11 


(11) 


Therefore,  if  the  induced  velocity  field  is  given  by  some  calculation  and  Ug^^ 
is  given  by  Eq.  (7)  or  Eq.  (8),  u^j^  can  be  solved  step  by  step  from  to 
Ugj.  Eq .  (6)  is  changed  to 


>’pi+i*  >'pi* 


<ii  .+ 

ni 


u  .+ 
ei 


u.  .+ 
11 


^i+l"^ 


ii+1 


(y 


ni+1 


(12) 


r  ...  can  be  solved  step  by  step  from  y  to  y  .  Eqs.  (11)  and 

pi+1  r  y  r  Pl  pn  ^ 

e  equations  to  be  solved.  However,  u^^^  Is  the  Induced  velocity  at 


Therefore,  y 
(12)  are  the^ 

JL 

the  position  y  .  and  y  .  is  an  unknown  quantity  till  Eq.  (12)  is  solved. 
Namely,  y  .  is  ^unknown ^when  Eq.  (11)  is  solved  at  first.  This  is  the  same 
about  y  ^when  Eq .  (8)  is  used.  Accordingly  at  first,  induced  velocity  and 

others  ire  calculated  by  assuming  'i  ^  After  solving  Eqs.  (11)  and  (12), 

the  calculated  y  .  is  used  for  ^recalculating  the  induced  velocity  and 
others.  The  sacSe'  operation  is  repeated  until  a  convergence  condition  is 
satisfied.  The  results  shown  in  section  3  are  those  for  five  iterations. 


Next,  the  method  of  taking  three-dimensional  correction  into  account  is 
shown.  Uhen  the  above  mentioned  equations  in  the  two-dimensional  case  are 
reviewed,  Eos.  (4),  (5)  and  (11)  are  entirely  the  same  as  those  of  Huang  et 
al.'s  method  for  the  axisymmetric  case.  This  is  the  natural  consequence, 
because  Eqs.  (^),  (5)  and  (11)  represent  the  condition  that  there  is  no  pres¬ 
sure  difference  through  vortex  layers  as  shown  in  Eq.  (10),  if  considering  the 
flow  field  with  discrete  vortex  layers.  Only  Eq .  (12)  which  determines  the 
location  of  a  vortex  layer  from  the  massflow  condition  is  different  from  that 
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of  the  axisymmetric  case.  Accordingly,  as  for  the  case  that  the  effect  of  the 
flow  contraction  in  the  z  direction  is  taken  into  account,  it  is  sufficient  to 
change  Eq.  (12)  to 


h„  (u  ,+  u  .  ,  , ) 

^  .  2nv  _ ni  ni+1 _ 

^pi+1  ^pi  n-  (u  .+  U..+  u  u. , , , 

^  ^  2p  ei  11  ei+1  11+1 


>'„i> 


(13) 


where,  h2  is  the  streamline  interval,  h^  and  h2jj  are  those  for  a  streamline 
passing  through  the  same  points  in  the  ^region  where  the  influence  of  the 
propeller  does  not  exist.  A  vortex  model  in  the  discrete  form  is  shown  in 
Fig.  4.  As  shown  in  the  figure,  the  velocity  distribution  on  a  line  of  z=con- 
stant  is  represented  by  the  distribution  of  m  and  the  movement  in  the  z 
direction  of  vortex  layers  due  to  propeller  suction  is  calculated  by  potential 
calculation.  At  this  time,  h2  /h2  corresponds  to  the  change  of  the  length  of 
a  vortex  filament.  Furthermore,  ^precaution  must  be  taken  for  the  fact  that 
the  z  coordinate  of  considering  velocity  distribution  with  propeller  is  dif¬ 
ferent  from  that  without  propeller.  Actually  h2  /^2n  considered  to  be  a 
function  of  y,  but  it  is  assumed  to  be  constant  bn  a  line  of  z=constant  for 
simplicity's  sake  in  this  discussion.  The  calculation  of  h2  is  carried 
out  in  the  uniform  flow.  The  uniform  flow  velocity  U  is  dete?'mined  by  thrust 
identity  from  the  calculated  thrust  and  the  calculated  propeller  open-water 
characteristics  similarly  to  the  self  propulsion  test.  That  is,  at  the  plane 
of  y  =  D/4, 


dz  _  1 

dx  U+u^ 


(14) 


is  traced  from  2D  ahead  of  the  propeller  and  the  movement  of  the  z  coordinate 
is  calculated,  where  D  is  the  diameter  of  a  propeller.  This  is  not  a  stream¬ 
line,  but  it  is  considered  that  the  flow  contraction  in  the  z  direction  can  be 
represented  to  some  degree  by  Eq.  (14),  which  is  consistent  with  the  fact  that 
the  present  method  is  a  simple  prediction  method. 

3.  Predicted  results  and  discussion 

3.1.  Flow  field  ahead  of  propeller 

To  examine  the  present  method,  the  total  velocity  distribution  with 
propeller  is  predicted  from  the  nominal  velocity  distribution  measured  ahead 
of  the  propeller,  and  compared  with  the  measured  total  velocity  distribution 
with  propeller.  The  experiment  was  carried  out  by  one  of  the  authors  [13], 
Principal  particulars  of  the  ship  model  and  the  propeller  are  shown  in  Table 
1.  The  comparison  is  carried  out  corresponding  to  the  ship  point  (ship  veloc¬ 
ity;  V  =  1.3m/sec,  number  of  revolution;  n=7.12  rps.).  Yanazaki's  method 
(infinitely  bladed  theory  [14])  is  used  for  the  propeller  calculation.  For 
the  calculation  of  flow  field  ahead  of  propeller,  induced  velocity  field  is 
calculated  by  using  uniform  flow  (l-w^)V  (w,p  is  the  effective  wake  fraction  by 
thrust  identity  method)  as  the  input.  Calculations  are  made  for  two  stations 
S.S.3/8  (x=-275mm)  and  S. S. 1 /4(x=-125nm) ,  where  velocity  measurements  were 
obtained  for  both  the  with-  the  without-propel ler  conditions. 
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First,  the  flow  contraction  in  the  z  direction  is  not  taken  into  account 
and  the  calculation  is  carried  out  by  use  of  Eqs.  (7),  (11)  and  (12).  The 
examples  of  the  calculated  results  for  station  S.S.  3/S  are  shown  in  Fig.  5. 
In  the  calculation,  the  dot-dash  lines  which  smoothly  connect  the  measured 
nominal  velocity  u^  plotted  in  circles  are  taken  as  the  input.  Then,  broken 
lines  of  effective  velocity  Ug  and  solid  lines  of  total  velocity  u  are  pre¬ 
dicted.  u  -  u  shown  with  dotted  lines  represents  the  velocity  increment  due 
to  propelfer  suction  and  the  difference  between  u  -  u  and  the  propeller 
induced  velocity  (u^  -  Ujj)  is  considered  to  show  thi  interaction  between  the 
propeller  and  the  stern  flow  field.  Since  the  induced  velocity  does  not 
change  a  lot,  it  is  considered  that  the  interaction  is  large  in  the  part  of 
large  Up-u^.  In  fig.  5,  a  fairly  good  agreement  between  measured  and  calcu¬ 
lated  total  velocity  distribution  is  obtained.  It  is  found  that  the  large 
effect  of  propeller  suction  is  observed  in  the  predicted  results  as  well  as  in 
the  measured  results  though  the  induced  velocity  is  very  small  in  the  station 
S.S. 3/8.  In  particular,  the  large  u  -  u  due  to  propeller  suction  is  pre¬ 
dicted  well  in  the  part  of  large  vilociEy  gradient  on  the  outside  of  the 
plateau-like  portion  of  the  velocity.  It  is  considered  that  the  large  veloc¬ 
ity  increment  due  to  propeller  which  cannot  be  predicted  well  by  an  ordinary 
boundary  layer  calculation  [13]  is  predicted  well  by  Eqs.  (7),  (11)  and 
(12).  Fig.  7  shows  how  the  effect  of  propeller  suction  appears  due  to  veloc¬ 
ity  distribution  in  y  direction.  Fig.  7(a)  shows  the  result  for  the  velocity 
distribution  around  a  fine  ship  like  the  Uigley  model.  It  is  found  similarly 
to  the  previous  paper  [7]  that  the  velocity  increment  becomes  larger  as  the 
point  gets  near  to  the  surface.  Fig.  (b)  shows  the  case  of  larger  velocity 
gradient  (vorticity  ^s  compared  with  (a).  The  velocity  increment  is 
larger  than  Fig  (a)  and  the  effective  velocity  distribution  becomes  close  to 
the  total  velocity  distribution.  However,  it  is  the  same  that  the  velocity 
increment  becomes  larger  as  the  point  gets  near  to  the  surface.  Fig.  (c) 
shows  the  results  for  the  example  of  the  velocity  distribution  around  a  full 
form  ship.  It  is  found  that  the  velocity  increment  distribution  is  consider¬ 
ably  different  from  the  case  of  (a)  and  (b).  The  examples  of  the  calculated 
results  for  the  station  S.Sl/4  are  shown  in  Fig.  7.  It  is  found  that  the  part 
of  large  velocity  gradient  on  the  outside  of  the  plateau-like  portion  of  the 
velocity  moves  too  inwardly  and  the  influence  of  propeller  suction  is  pre¬ 
dicted  to  be  too  strong.  However,  it  seems  that  the  qualitative  tendency  of 
the  velocity  increment  distribution  agrees  well  with  experimental  results.  In 
order  to  show  clearly  this  feature,  the  comparison  between  measured  and  pre¬ 
dicted  velocity  increment  distribution  is  shown  in  Fig.  8  and  Fig.  9.  As 
shown  in  Fig.  8,  the  predicted  result  agrees  well  witli  measured  ones  such  as 
in  Fig.  5  in  S.S.  3/8,  where  the  induced  velocity  is  small.  As  for  the  result 
in  S.S. 1/4  shown  in  Fig.  9,  it  seems  that  the  predicted  position  of  the  region 
where  the  velocity  increment  is  remarkable  agrees  well  with  measured  one, 
though  the  magnitude  of  velocity  increment  is  larger  by  about  0.1  than  the 
measured  one.  Therefore,  it  is  considered  that  even  the  simple  method  of 
using  Eqs.  (7),  (11)  and  (12)  can  be  used  for  judging  whether  the  ratio  (1- 
W-^)  /  ( l"On)nieari  large  or  small  owing  to  the  characteristic  of  nominal  veloc¬ 
ity  distribution  ((^“'^n^mean  “  ^^n^mean’  volume  mean  of  nominal  velocity 
in  the  propeller  disk,  W^;  effective  wake  fraction). 

Next,  Fig.  10  shows  the  results  for  the  station  S.S. 1/4  predicted  by  the 
present  method  using  Eqs.  (8),  (11)  and  (12),  in  which  the  effect  of  the  flow 
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contraction  in  z  direction  and  the  change  in  potential  wake  for  the  outside 
edge  condition  are  taken  into  account.  The  potential  velocity  is  calculated 
by  the  Hess-Snith  method.  As  shown  in  Fig.  10,  a  fairly  good  agreement 
between  measured  and  predicted  total-velocity  distribution  near  the  propeller 
axis  is  obtained  as  a  whole,  though  disagreement  is  somewhat  large  at  some 
places  near  the  plateau-like  portion  of  the  velocity.  The  velocity  contour 
curves  in  S.S.1/4  are  shown  in  Fig.  11.  On  the  right  side  of  the  figure,  the 
comparison  between  the  measured  and  predicted  total  velocity  distribution  is 
shown.  On  the  left  side,  the  comparison  between  the  measured  nominal  velocity 
distribution  and  the  predicted  effective  velocity  distribution  is  shown.  From 
the  figure  of  the  total  velocity,  the  measured  and  the  predicted  results  agree 
fairly  well  with  each  other  as  a  whole,  although  the  predicted  velocity  is 
somewhat  higher  than  the  measured  one  near  the  upper  and  lower  ends  of  the 
propeller  disk.  That  the  effect  of  the  flow  contraction  in  the  z  direction  is 
calculated  in  the  uniform  riow  of  (l-w^)V  for  simplicity's  sake  and  the  irro- 
tational  velocity  component  in  the  boundary  layer  is  not  taken  into  account  in 
the  similar  form  to  that  of  the  induced  velocity  when  deriving  Eqs.  (l)-(4)  is 
considered  to  be  the  reason  for  the  discrepancy.  It  also  seems  that  the  neck¬ 
like  portion  of  the  nominal  velocity  distribution  owing  to  the  cross  flov; 
appears  too  strongly  in  the  predicted  total  velocity  distribution.  Other 
methods  seem  to  be  necessary  for  the  inproveiaent  in  this  part  as  the  effect  of 
the  cross  flow  is  not  taken  into  account.  However,  it  is  considered  that  the 
effect  of  propeller  suction  for  the  flow  field  around  a  full  form  ship  is 
predicted  fairly  well  as  a  whole  by  present  method,  if  considering  that  the 
results  in  Fig.  10  and  Fig.  11  are  predicted  by  the  method  assuming  that  the 
dominant  effect  of  propeller  suction  on  the  stern  flow  field  can  be  treated 
two-dimensional ly.  From  the  left  figure  of  Fig.  11  and  Fig.  10,  the  differ¬ 
ence  between  the  nominal  and  the  calculated  effective  velocity  distribution  is 
investigated.  It  is  observed  that  the  contour  curves  of  the  effective  veloc¬ 
ity  distribution  are  displaced  toward  the  center  plane  as  compared  with  those 
of  the  nominal  velocity  distribution  as  a  whole  and  the  velocity  is  increased 
from  the  nominal  one  to  the  effective  one,  but  not  uniformly.  Namely,  there 
are  the  regions  where  the  nominal  and  Ihe  effective  velocity  are  hardly  dif¬ 
ferent  and  the  regions  where  the  difference  between  them  is  large.  Therefore, 
it  .3  found  that  the  effective  velocity  distribution  by  present  method  is 
considerably  different  from  the  distribution  determined  from  the  nominal 
velocity  distribution  being  multiplied  by  ^  ^  ^~^^n^mean»  which  has  been 
used  frequently  for  propeller  calculations. 


Next,  in  order  to  investigate  whether  the  information  can  be  obtained  or 
not  about  the  change  of  velocity  in  a  cross  section  (Vg-V^,  the  dis¬ 

tortion  of  the  longitudinal  vorticity  distribution  due  to  propeller  suction  is 
predicted.  The  longitudinal  vorticity  with  propeller  is  predicted  from  the 
measured  one  without  propeller  by  the  method  taking  into  account  only  the 
decrease  of  the  cross  sectional  area  and  the  change  of  position  of  stream 
tubes,  which  are  obtained  as  the  results  of  prediction  of  the  propeller  effect 
for  the  velocity  in  x  direction  by  present  method.  Namely,  if  the  longitu¬ 
dinal  vorticity  without  propeller  is  known,  the  longitudinal  vorticity 

-1)^  v;ith  propeller  is  predicted  by 
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where  Che  pareachesized  terms  show  the  position.  The  comparison  of  Che  pre¬ 
dicted  and  measured  longitudinal  vorticity  is  shown  in  Fig.  12.  The  predicted 
results  agree  well  with  measured  ones  with  propeller  as  a  whole,  except  near 
the  surface.  Eq.  (15)  is  able  to  predict  the  increase  of  the  peak  value,  the 
raoveraenc  of  the  peak  position  and  the  tendency  that  the  width  of  vorticity 
distribution  becomes  narrow  when  Che  propeller  is  operating.  But,  it  is  found 
that  the  distribution  is  contracted  somewhat  excessively  toward  the  surface  on 
the  outside  of  the  peak.  The  contour  curves  of  the  predicted  vorticity  dis¬ 
tribution  with  propeller  is  compared  with  the  measured  vorticity  distributions 
with  and  without  propeller  in  Fig.  13.  It  also  seems  that  the  prediction 
method  is  able  to  predict  the  peak  value  and  position  to  some  extent.  There¬ 
fore,  it  is  considered  that  the  information  regarding  V^-  can  be 

obtained  to  some  degree  by  present  method  from  the  difference  of  the  velocity 

induced  by  m  and  oi 

xn  xp 

3.2  Effective  velocity  distribution  on  propeller  plane 

The  effective  velocity  distribution  on  the  propeller  plane  is  predicted 
from  the  measured  nominal  velocity  distribution.  In  this  calculation,  iter¬ 
ative  calculation  is  needed  to  calculate  the  strength  of  the  bound  vortex  in 
the  effective  velocity.  The  induced  velocity  is  calculated  at  the  first 
iteration  by  using  the  nominal  velocity  distribution  as  the  input  and  there¬ 
after  the  predicted  effective  velocity  distribution  is  used  as  the  input  to 
recalculate  the  induced  velocity.  The  sane  operation  is  repeated  until  a 
convergence  condition  is  satisfied.  In  this  paper,  four  iterations  were 
needed.  Since  y^  is  taken  as  the  center  line  in  Eq.  (9)  and  the  calculation 
is  carried  out  only  for  starboard  side  of  the  ship,  the  induced  velocity  at 
the  position  D/4  ahead  of  the  propeller  is  used.  The  comparison  between  the 
measured  nominal  and  the  predicted  effective  velocity  distribution  is  shown  in 
Fig.  14.  The  right  side  is  the  measured  nominal  velocity  distribution  used  as 
the  input.  (Measurement  was  carried  out  for  the  sane  ship  and  sane  advance 
speed  by  Tanaka  et  al.  [15],  but  the  data  analysis  method  for  the  five-hole 
pitot  tube  is  somewhat  different  from  that  of  Kasahara  [13]).  The  left  side 
is  the  predicted  effective  velocity  distribution.  It  is  found  similarly  to 
3.1  that  the  contour  curves  of  the  effective  velocity  distribution  are  dis¬ 
placed  toward  the  center  plane  as  compared  with  those  of  the  nominal  velocity 
distribution  as  a  whole  and  the  velocity  is  increased  from  the  nominal  one  to 
the  effective  one,  but  not  uniformly.  The  velocity  increment  from  the  nominal 
to  the  effective  velocity  is  large  particularly  on  the  outside  of  the  plateau¬ 
like  portion  of  the  velocity.  When  focusing  on  the  propeller  disk,  the  part 
of  large  velocity  gradient  is  hardly  in  it  in  the  case  of  the  nominal  velocity 
distribution,  while  high  velocity  region  enters  into  it  in  the  case  of  the 
effective  velocity  distribution.  For  reference,  the  volume  mean  of  effective 
velocity  is  compared  with  that  of  the  nominal  velocity  and  the  mean  inflow 
velocity  from  the  self-propulsion  test  on  the  basis  of  the  thrust-identity 
method  in  Table  2.  The  calculated  thrust  by  using  respective  velocity  distri¬ 
butions  as  the  input  are  also  shown  in  Table  1.  It  is  found  that  the  present 
method  is  able  to  predict  the  tendency  in  the  change  of  the  flow  field  due  to 
propeller.  But  the  velocity  increment  from  the  nominal  velocity  to  effective 
one  is  considered  to  be  predicted  somewhat  excessively.  Since  the  predicted 
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total  velocity  is  sonewhat  higher  than  the  measured  one  at  some  places  in  the 
results  shown  in  3.1,  the  reason  for  the  slightly  excessive  prediction  in  this 
section  is  considered  to  be  similar  to  that  discussed  in  3.1.  However,  as  the 
nominal  velocity  distribution  was  measured  very  much  before  the  self  propul¬ 
sion  test,  further  discussion  is  withheld.  It  is  considered  that  the  change 
from  the  nominal  velocity  to  the  effective  one  is  predicted  to  some  degree 
from  the  results  in  3.1  and  in  Table  2.  Therefore,  in  order  to  investigate 
the  characteristics  of  the  effective  velocity  distribution,  the  distributions 
of  the  circumferential  mean  nominal  and  effective  velocity  in  the  radial 
direction  are  shown  in  Fig.  15  and  the  distributions  of  nominal  and  effective 
velocity  in  circumferential  direction  at  three  r/R  are  shown  in  Fig.  16  (r; 
distance  from  the  propeller  axis,  R;  propeller  radius).  From  Fig.  15,  it  is 
found  that  the  circumferential  mean  of  the  nominal  velocity  varies  gently  in 
the  radial  direction,  and  the  circumferential  mean  of  the  effective  velocity 
is  increased  as  a  whole  due  to  propeller  suction  from  that  of  the  nominal 
velocity,  especially  the  velocity  increment  is  large  near  the  propeller  hub 
and  the  blade  tip.  In  particular,  near  the  blade  tip,  the  velocity  increment 
is  large  by  pulling  the  high  velocity  region  into  the  propeller  disk,  as  men¬ 
tioned  about  Fig.  14.  Therefore,  the  predicted  velocity  increment  distribu¬ 
tion  due  to  propeller  for  the  full  form  ship  is  different  from  that  for  a  body 
of  revolution  [4]  and  that  for  a  simple  Uigley  model  [7].  Though  the  propel¬ 
ler  is  relatively  small  in  the  ship  model  used  in  this  paper  and  the  high 
velocity  region  of  nominal  velocity  distribution  on  the  outside  of  the  pla¬ 
teau-like  portion  hardly  enters  into  the  propeller  disk,  it  is  considered  that 
the  tendency  of  the  large  velocity  increment  being  observed  on  the  outside  of 
the  plateau-like  portion  is  the  same  for  the  full  form  ship  having  a  larger 
propeller.  Therefore,  it  is  considered  that  the  tendency  of  change  from  the 
nominal  velocity  distribution  to  the  effective  velocity  distribution  for  a 
full  form  ship  is  considerably  different  from  that  for  a  body  of  revolution  or 
a  fine  ship.  Also  at  r/R  =  0.9813  near  the  blade  tip  in  Fig.  16,  the  effec¬ 
tive  velocity  is  remarkably  large  as  compared  with  the  nominal  velocity  in  the 
region  60  <  9  <  150  and  the  unevenness  of  the  effective  velocity  distribution 
in  circumferential  direction  is  increased  from  that  of  the  nominal  velocity. 
This  phenomenon  for  the  full  form  ship  is  different  from  that  for  a  simple 
Wigley  model,  in  which  the  velocity  increment  from  the  nominal  velocity  to  the 
effective  one  becomes  larger  as  the  point  gets  near  to  the  surface. 

4.  Conclusion 

The  distortion  of  stern  flow  field  around  a  full  forra  ship  due  to  propel¬ 
ler  suction  is  investigated.  The  effect  of  propeller  suction  is  predicted  by 
a  modified  form  of  Huang  et  al. 's  method  under  the  assumption  that  the  effect 
of  the  change  of  vorticity,  vj'nose  axis  is  vertical,  is  dominantly  responsible 
for  the  distortion  of  velocity  distribution  due  to  propeller  suction.  The 
conclusions  are  as  follows. 

(1)  The  simple  prediction  method  in  this  paper  is  able  to  properly  repre¬ 
sent  the  effects  due  to  propeller  suction  observed  in  the  measured  results. 

(2)  The  change  of  longitudinal  vorticity  distribution  can  be  predicted  by 
taking  into  account  the  decrease  of  the  cross-sectional  area  and  the  change  of 
position  of  stream  tubes,  which  appear  in  present  prediction. 

(3)  The  effective  velocity  around  a  full  form  ship  is  remarkably  increased 
in  the  region  of  large  velocity  gradients  on  the  outside  of  the  plateau-like 
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portion  of  the  velocity  from  the  nominal  velocity.  From  this  fact,  the 
unevenness  of  the  effective  velocity  distribution  in  circumferential  direction 
is  increased  from  that  of  the  nominal  velocity  neur  the  blade  tip  around  a 
full  form  ship. 
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Table.  1  Principal  Dimension  of  model 


Lpp 

1  2  m 

8 

2.182  m 

Q- 

d 

0.  726  m 

<7) 

Cb 

0.837 

Rn 

1.64x10^ 

0 

0.301  m 

LU 

Pitch  ratio 

0.6155 

O. 

o 

EAR 

0.6050 

Boss  ratio 

0.1606 

Table  2  Comparison  of  thrusts  and  mean  velocity  (n=7.t2  rps) 


- 1 

input  velocity 
distribution 

nominal  velocity 

Un 

effective  velocity 

u» 

uniform  flow 
(l-wi)V 

thrust 

8.67 

7.77 

8.05 

lUnlv 

(u«lv 

l-wt 

lUn)o.7R 

mean  velocity 

.393  (.382)* 

.494 

.461 

.388 

mm 

[  Jv  denotes  volume  meon  of  velocity  distribution 
(  lo.7R  denotes  mean  velocity  of  0.7R 
•  is  the  (unjv  calculated  at  NKK 


wi  denotes  effective  woke  froction  by  thrust  identity 
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Fig'  n  Coapftriaon  of  velocity  diatributioo 
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